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Abstract 

In this paper we study the time-dependent Ginzburg-Landau equations on a smooth, 
bounded domain 1] C R 2 , subject to an electrical current applied on the boundary. The 
dynamics with an applied current are non-dissipative, but via the identification of a special 
structure in an interaction energy, we are able to derive a precise upper bound for the energy 
growth. We then turn to the study of the dynamics of the vortices of the solutions in the limit 
e — > 0. We first consider the original time scale, in which the vortices do not move and the 
solutions undergo a "phase relaxation." Then we study an accelerated time scale in which the 
vortices move according to a derived dynamical law. In the dynamical law, we identify a novel 
Lorentz force term induced by the applied boundary current. 



o\ : 1 Introduction 

o ■ 

1.1 Formulation of the equations and boundary conditions 



The Ginzburg-Landau free energy functional on a domain Q C M 2 is denned for a function u : Q — > 



C, a vector field A : Q — > IR 2 , and a real parameter e > by 



\ 2 \2 



F £ (u, *) = \J 1 VaM ' 2 + |CUrl + 2e" ' ' (L1) 

where Vau '■= Vw — iAu is the covariant gradient of u. Ginzburg and Landau introduced their 
eponymous functional in 1950 [IE] as a free energy in a phenomenological model of superconductivity. 
In this setting Q is thought of as a two-dimensional cross section of a sample of superconducting 
material; we shall assume Q is smooth and bounded. The function u is known as the "order 
parameter," and models the relative density and phase of superconducting electrons, with \u\ ~ 1 
indicating the superconducting state and |u| ~ indicating the normal state. The vector field 
A is the magnetic vector potential and h := curl A is the induced magnetic field strength in Q. 
The vector he^, which is orthogonal to Q, is the induced magnetic field. In the model, e > is 
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a parameter depending on the material comprising the superconductor. We are interested in the 
regime £ < 1, which corresponds to so-called "extreme type-IP superconductors. 

The equations modeling the dynamics of superconductors, derived by Gor'kov and Eliashberg 
in 1968 [T7], are the covariant heat flow for the Ginzburg-Landau functional: 

d t u + i$u = A A u + ±(l - \u\ 2 ) in x 1+ 

a(d t A + V$) = V ± h + (iu, V A u) in x R+. ^ ' ' 

In these equations we have written A A u '■= (div — iA-)V au for the covariant Laplacian and V x h := 
(—d2h,dxh) for the perpendicular gradient of the induced magnetic field. The function $ : Q — y R 
is the electric potential, and E := — (d t A + V$) is the induced electric field. The constant a > is 
the conductivity of the superconducting material, and by Ohm's law the quantity oE is the current 
of normal (i.e. non-superconducting) electrons in Q. The vector field (iu, V^ii) is known as the 
supercurrent, and represents the current of electrons in the superconducting state. Here we have 
employed the notation (a, b) := 9ft(a)9?(5) + Q(a)Q(b) for a, b G C to mean the inner-product with 
C identified with R 2 , and (a, X) for X G C 2 to mean the vector in M 2 with components (a, Xi) and 
(a,X 2 ). 

The equations ( II. 2p give rise to the Maxwell equations in Q. The second equation in (II. 2p is 
Ampere's law with the total current given as the sum of the normal and super currents. Faraday's 
law of induction is seen in the relation curl-E = — curl(<9jv4 + V$) = —dfh. Taking the inner- 
product (iu, ■) with the first equation in (jl.2p and taking the divergence of the second, we find 
div E = —(iu,d<$,u) /a, where we have written d^u = d t u + iQu. Then —(iu,d$u)/a is the charge 
in Q, and Gauss's law holds. The remaining Maxwell equation, Gauss's law of magnetism, follows 
since div(/ie 3 ) = d 3 h = 0. 

The evolution of the superconductor is coupled to the electromagnetic fields on the exterior of 
the domain, ST = R 2 \fi, by assuming that the electromagnetic fields satisfy Maxwell's equations 
everywhere [T71 [HI [13]. In the absence of surface charges, this gives rise (cf. [18]) to boundary 
jump conditions coupling the external electric and magnetic fields, E ex and H ex , to the fields in Q. 
These are 

E ■ v = E ex ■ v and h = H ex , (1-3) 

where v is the outward unit normal on <9f2. We will assume that the exterior fields satisfy the static 
Maxwell equations, which in particular requires 

V ± H ex = -I ex , (1.4) 

where I ex : Q c — > IR 2 is the smooth, time-independent current of normal electrons on the exterior. 
If Q c is a conductor, we will assume it has the same conductivity as Q so that I ex = uE ex . If Q c 
is not a conductor, then I ex = and we take H ex > 0; E ex need not vanish, but it is frequently 
assumed to when I ex does. 

It is mathematically convenient to recast the boundary conditions for E ■ v in (I1.3P in terms of 
a condition for V ' au ■ v. Since V L h ■ v = Vh ■ r = —d T h with r the unit tangent, we may plug 
h = H ex and (11.41) into the second equation in (II .2p to see 

(iu, V au ■ v) = I ex ■ v — aE ■ v. (1.5) 

When oE ex = I ex (with I ex possibly 0) we may then take V^m- v = 0, which implies the appropriate 
boundary condition, E ex ■ v = E ■ v. 
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We will actually study a generalization of these boundary conditions, which for clarity we record 
now along with the evolution equations: 

d t u + i<$>u = A A u + ^(1 - \u\ 2 ) in n x 1+ 

d t A + V$ = V x h + (iu, V A u) in O x 1+ 

Vau ■ v = iuJ ex ■ v on dVl x 1R + (1.6) 

h = H ex on dn x R+ 
(u, A$)|t=o = (u ,A ,%). 

Here H ex is again given by (jl.4p , but now we take J ex : dQ — > M 2 to be any smooth vector field. To 
reduce notational clutter, we have assumed that a = 1, but all of our results may be modified to 
handle any fixed a > 0. The introduction of J ex is justified in three ways. First, from a mathematical 
point of view, the case J ex ^ is just a generalization of the case J ex = 0. The methods we develop 
to handle — V" L i? ea; = I ex ^ also handle J ex ^ 0, which justifies referring to J ex as a sort of current. 
Second, the actual physical jump condition across dfl is [18J 



E ■ u = E ex ■ v + q, (1.7) 

where q is the surface charge accumulated on dfl. Plugging the generalized condition V^m • v = 
iuJ ex ■ v into ( II. 5p when I ex = E ex yields 

E ■ v = E ex ■ v - \uf J ex ■ v, (1.8) 

which shows that J ex ■ v behaves as a sort of surface charge. Similarly, when I ex = but E ex ^ 0, 
which corresponds to an external voltage, ( II. 5p gives 

E ■ v = — \u\ 2 J ex ■ v (1-9) 

so that J p x ' v can also behave as an external voltage would. This suggests that the generalized 
boundary condition can be used as an approximate model of surface charge or external voltage. 
Third, the currents J ex and I ex are independent, which gives a mechanism for inducing different 
scales of current forcing. The reader solely interested in the standard choice of boundary conditions 
h = H ex , Vau ■ v — may simply take J ex = in all of our analysis. 

The novelty in our analysis is in the presence of the applied currents I ex and J ex . Numerous 
authors [3j |23j (2TJ [191 El El 13 1211 [29] have studied the non-magnetic analog of the equations (II. 6p 
for which an applied current is impossible. In the magnetic case (II. 6p has been studied rigorously 
with I ex = J ex = and H ex a constant in [311 [26]. Several numerical and formal asymptotic results 
are available [T21 [TBJ El] when I ex ^ and J ex = 0, a stability analysis of the normal state (u = 0) 
with applied current was performed in [2], and a 1-D model of a superconducting wire with current 
was studied in [21], but we are aware of no rigorous results in the 2-D magnetic model with applied 
current or with a surface charge. 



1.2 Definitions and terminology 

We will now record several definitions and bits of terminology that will be used throughout the 
paper. For a more thorough exposition of these quantities and of the magnetic Ginzburg-Landau 
model in general, we refer to the book [27J and the references therein. 

The objects of interest in the study of ( II. 6p are the zeroes of the complex-valued function u; 
these are known as vortices. Each vortex carries an integer topological charge called its degree, 
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which is defined as the winding number of the map uj \u\ on any simple closed curve around the 
zero. The energy density (the integrand of F e ) concentrates around the vortices, with the i th vortex 
contributing an amount of energy of the order n \di\ |loge| to F £ , where di is the degree of the vortex. 

The energy F e possesses a gauge invariance under the pointwise action of the group U(l): 
F e (u,A) = F £ (ue l ^,A + V£) for any sufficiently smooth £ '■ — > K. This gauge invariance carries 
over to solutions of the equations (11. 6p as well in the sense that if (u, A, $) are solutions, then so 
are (ue^, A + V£, $ - d t £) for £ : x R+ -> R. We refer to the change u H- ue^, A \-t A + V£, 
$ i— t- $ — <9 t £ as a gauge change. For solutions to (jl.fip to be unique, we must eliminate the gauge 
invariance by "fixing a gauge." In this paper, we work exclusively in what we call the $ = / gauge 
(see Lemma [2.41) . 

For any e > and any choice of J ex , H ex smooth, the system (I1.6P is well-posed for all time in 
any fixed gauge, and the solutions are smooth. This may be established through a modification of 
standard results [TD1 [TTJ . For the sake of completeness, in Appendix |A] we make a few brief remarks 
on the necessary changes and the a priori estimates available for solutions. 

The vortices of a configuration (u,A) are best described through the "vorticity," fi(u, A), a 
gauge-invariant version of the Jacobian determinant of u: 

fi(u, A) := curl(ra, V au) + curl A. (1-10) 

It is now well-known in the literature (cf. [20] for first results but without magnetic field, and 
[271 S3] for results with magnetic field) that fi(u, A) ~ 27r^<ij<5 ai , where G Q, di G Z are the 
location and degree of the i th vortex. Here ~ means close in various norms: (C 0,a (Q))* , for instance. 

As in the non- magnetic case [1], two mechanisms contribute to the energy F e . In the case we 
are interested in, with n vortices of degree di = ±1, the result is that (roughly speaking) 

F E (u, A) = n(7i \loge\ + 7 ) + W d {a) + 0(1). (1.11) 

The first term on the right is the self-energy of the vortices, which is itself composed of a divergent 
term 7m|loge| and a finite term nj, where 7 is a known constant related to the structure of a 
vortex. The second term is the inter-vortex interaction energy, or magnetic renormalized energy. 
For a collection of points a = (ai, . . . , a n ) G fl n and degrees d = (di, . . . , d n ) G Z n , the function 
Wd(a) is defined by 

Wd(a) = —ix didj log \di — dj \ + 7r didjSa(ai, aj). (1-12) 

Here Sq G C l (Q x Q) is the regularization of the Green's function for the London equation on Q, 
i.e. 

S a (x,y) = G(x,y)+log\x-y\, (1.13) 

where 

-A x G(x, y) + G(x, y) = 2n5 y in Q 
G(x,y) = for x EdQ,y EQ. 

See |27J for further discussion of the magnetic renormalized energy and for proof of (11.111) . 

Throughout the paper we will use the notation o(l) to refer to any quantity that vanishes as 
£ — > 0. Similarly, 0(1) refers to a quantity that stays bounded. For two quantities a £ , b £ , we will 
employ the notation a £ b £ to mean that a £ /b £ = o(l), and we write a £ x b £ if a £ /b £ = 0(1) and 
b £ /a £ = 0(1). We will also employ the standard convention of using the letter C to denote a generic 
positive constant that may change from from line to line. 
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1.3 Known results and expectations 



The main interest in studying (11. 6ft is to derive the dynamics of the vortices associated to sequences 
of solutions (u £ , A £ ,$> £ ) in the limit as e — > 0. In the non-magnetic case, this was accomplished 
under various assumptions and for varying lengths of time [31 [23j [211 HI El El [28], [29]. The 
magnetic case with I ex = J ex = and H ex = h ex (e) a constant depending on e was studied in [31] 
for h ex (e) fixed and in J2H] for h ex (e) = (3 \loge\ with (3 > 0. In this setting the equations (jl.fip 
constitute the L 2 gradient flow of the full Ginzburg-Landau energy, 

G e (u, A) = ±J \V A u\ 2 + \h- h ex (e)\ 2 + (1 ~^ |2)2 . (1.15) 

This leads to energy dissipation, G e (u e (t), A e (t)) < G e (u s (0),A E (0)) for all t > 0, which provides 
the technical advantage of a priori control of the energy and the number of vortices. 

The general scheme found in these papers is that the vortices do not move until after an amount 
of time of order A £ = |loge| /h ex (e). When h ex (e) is fixed, this means that in the limit e — > the 
vortices cannot move at all since A e — > oo. In this case (u £ ,A e ,§ e ) — > (u*, A*, $*) in some sense, 
and it is possible to pass to the limit in (11. 6p to derive the dynamics for (u* , A* , <&*) . To see vortex 
motion, the solutions are accelerated in time at scale X e according to 

u E (x, t) ^ u e (x, X E t), A e (x, t) ^ A e (x, X £ t), $ e (x, t) ^ $ e (x, X s t). (1.16) 

In this scaling, the vortices do move in the limit, and their dynamics are governed by 

cn(t) = —V ai W d (a(t)) - 2d i h ex VH (a i (t)), (1.17) 

7T 

where aj(t) G Q is the location of the i th vortex, di G Z is its degree, V 0j Wd is the derivative of 
magnetic renormalized energy ( 11. 12ft with respect to a, G Q, and Vi^o is a purely magnetic forcing 
term with Hq the solution to the London equation 

-A# + #o = in ft 
i7 = 1 on 9ft. 

As such, the limiting dynamics are a gradient flow of the energy 

n 

W d , hex (a) := ^(a) +2irh ex ^^(fli), (1.19) 

i=l 

the latter term of which is the interaction energy between the vortices and the applied magnetic 
field. 

In [26], the choice h ex = (3 |loge| implies A £ = 1//3, which allows the vortices to move in the 
original time scale. The interpretation of this is that the interaction energy between the vortices 
and the magnetic field is sufficiently strong to induce vortex motion in the original time scale. 
The resulting motion corresponds to the gradient flow of this interaction energy, and is written 
dj = —2(3diVHQ(ai). Note that in this case the vortices do not interact with each other in the sense 
that the motion of the point aj does not depend on the points aj for j ^ i. 

In the case J ex ^ or I ex ^ the gradient flow structure of the equations (jl.6p breaks down. 
Energy does not dissipate, and we can no longer expect the limiting dynamics of the vortices to be a 
gradient flow. This creates serious difficulties in applying the standard Ginzburg-Landau toolboxes, 
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which rely crucially on precise knowledge of the energy. In fact, the applied boundary current is 
expected to introduce two novel features to the dynamics, both of which have been observed in 
the numerical simulations of [121 EH EH]- First, the applied current generates an electric field in 
Q, and the vortices feel a Lorentz force perpendicular to this field. Second, and more drastic, a 
sufficiently strong applied current is expected to create and destroy vortices near the boundary, 
thereby injecting or removing large amounts of energy from the system. 

1.4 Summary of main results 

In this paper we analyze sequences of solutions (u £ , A e , $ £ ) to ( 11. 6ft as e — > in both the original 
and accelerated time scales. Our aim is to show that the applied boundary currents I ex and J ex 
induce Lorentz forcing terms in the limiting vortex dynamics for the accelerated time scale, and to 
identify the structure of the Lorentz forces. We make the structural assumption that 

Jex = jexJy Hex = h ex H, and I ex = h ex I (1-20) 

for field strengths j ex = j ex (e) > and h ex = h ex {e) > 0, J : dfl — > M. 2 a smooth, fixed vector field, 
and H : dQ — > R the smooth trace onto dfl of the solution to the static exterior Maxwell equations 
V X H = —I. When I = we will assume that H = 1, corresponding to a uniform applied magnetic 
field. For notational convenience we will not write the e dependence for the parameters j ex or h ex . 
We shall consider four distinct regimes for the parameters: 
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(1.21) 



The first regime handles the choice of any j ex and h ex fixed with respect to e since we may simply 
rescale J and / to set j ex = h ex = 1. In the second two cases, at least one parameter blows up, 
but one dominates the other. In the fourth case both blow up but are of the same order. The 
upper bound by [logel^ 9 is purely technical, being required in the proof of the dynamical law in 
the accelerated time scale. We define the dominant field strength via 

k ex := max{j ex , h ex }- (1-22) 

The main thrust of the paper is to deal with the complications caused by the dynamics no longer 
being dissipative or even conservative (as with the Schrodinger flows associated to F e ). We show 
that, while the energy does not necessarily decrease, it cannot increase too quickly. This allows us 
to identify a time scale depending on k ex in which the number of vortices remains constant, but the 
vortices move in Q, exhibiting the additional Lorentz force drifts due to the applied currents. The 
control of the energy growth is far from trivial; indeed, we see that standard tricks (e.g. Gronwall) 
are insufficient for getting precise estimates of the energy. It is only via the identification of some 
very special structure in the current forcing terms that we are able to get the delicate estimates 
required. 

To understand how the free energy of solutions evolves in time, we introduce a splitting of 
the solutions into a topological (i.e. generated by vortices) component and an applied current 
component. This is similar to a technique employed in [SI El E] for the non-gauged case. To 
motivate the splitting, we take the curl of the second equation in ( II. 6p to see that 

d t h £ - Ah £ + h e = fx(u E , A £ ) in Q, ^ 
h e = h ex H on dVt. 
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We can then split h £ according to h £ = h ex h + h' £ where 

-Ah + h = mfl 
ho = H on d£l. 

Then ho is the static contribution of the applied magnetic field, and h' e , which satisfies h' e = on 
dfl, is the dynamic part of the induced magnetic field generated by the vortices. Notice that 

h' e = h £ — h ex h = curlv4 £ — h ex Ah = curlv4 £ — h ex curl V ± h = curl(v4 £ — h^V^ho). (1-25) 

This suggests defining the modified vector potential B £ := A e — /i e;E V~ L /io, which we expect to be 
the part of the vector potential generated by vortices. Then on dfl 

V Bs « £ ■ v = Va,u £ ■ v + iu £ h ex SJ^ho ■ v = iu £ J ex ■ v + iu £ V^E ex ■ v = iu £ (J ex - I ex ) ■ u, (1.26) 

which shows that J ex acts as a sort of current. It is then useful to modify u £ in such a way to turn 
the inhomogeneous Neumann boundary condition into a homogeneous one. To do so we define f\ 
and f as the solutions to 

' -A/i + A = -Afo + fo = inQ 

\7fx-v = J-v on 90 (1.27) 

V/o ■ v = I ■ v on dQ 

and write 

/ := JeJi - h. rU (1.28) 
Then the modified order parameter v £ := u £ e~ l f satisfies the homogeneous boundary condition 

Vb £ u £ ■ v = e~ if V Be u £ ■ v - iv £ Vf ■ v = iv £ (J ex - I ex - V/) • v = 0. (1.29) 

It will be convenient to introduce the forcing vector field 

Z £ ■= JsrV/i - h ex Vf - h ex V ± h . (1.30) 

In studying v £ ,B £ , we are led by the equations (11.61) to consider the evolution of a modification 
of the standard Ginzburg-Landau free energy F £ : 

F £ (v £ ,B £ ):=F £ (v £ ,B £ ) + l [ \v £ \ 2 \Z £ \ 2 . (1.31) 

The study of F £ (v £ , B £ ) points to a natural choice of gauge: one in which <3> £ = /, which we refer to 
as the $ = / gauge. 

Many of the arguments in this paper rely crucially on the initial data satisfying a well-preparedness 
condition on the energy. In particular, for a sequence (v £ , B £ ) and a constant Co > 0, we say that 
initial data (v E (0), B e (0)) are well-prepared at order C if /i(t> £ (0), B £ (0)) — > 2-n Y^=i di<5 ai with 
di = ±1 and 

F £ (v £ (0),B £ (0))<nn\\og6\ + W d (a)+ ni + l [ \Z £ \ 2 + C k eX} (1.32) 

1 Jn 

where Wd(a) is the renormalized energy defined by fll.12)) . Z £ is given by fll.30p . k ex is defined by 
( 11.221) . and 7 is a fixed constant (see Lemma \B. II) . We note that by adapting results in the literature 
( [27] for example) we may construct initial data satisfying these hypotheses. 
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The utility of studying v e ,B £ rather than u e ,A e lies in a novel observation on the structure of 
a term arising in the equation for the evolution of the modified energy. We find that if (u £ , A £ , $ e ) 
solve (jl.6p in the $ = / gauge, then for v £ = u £ e~ l f , B £ = A £ — h^V^ho, 

d t F £ (v £ ,B £ ) + [ \d t v £ \ 2 + \d t B £ \ 2 = [ V £ -Z £ , (1.33) 
Jn Jn 

where V £ is the "velocity component" of the full space-time Jacobian associated to (v e , B £ ) (see 
Section [2]). The extra structure of the interaction term, V £ ■ Z £ , is the key to controlling the growth 
of the modified energy because of estimates for V £ proved in [25] (recorded here in Proposition 12. 14j) . 
Using these estimates, we can prove that if the initial data is well-prepared at order Co, then in 
an amount of time of order A e := ^? g ^ , the modified energy F £ (v £ ,B £ ) can increase at most by an 
amount 2C§k ex . 

Theorem 1 (proved later as Theorem I2.16p . Suppose that the initial data (v £ (0), B £ (0)) are well- 
prepared at order Cq, as defined by (11.321) . Let \ £ = |loge| /k ex . Then there exists a constantT > 
so that, as e — > 0, 

F £ (v £ ,B £ )(t) < F £ (v £ ,B £ )(0) + 2C k ex for all t e [0,T A £ ], and (1.34) 

/■ToA £ r- 

/ / \d t v £ \ 2 + \d t B £ \ 2 <2C k ex . (1.35) 
Jo Jn 

When j ex = h ex = 1 we can derive the limiting dynamics in the original time scale for the pair 
(v £ , B £ ) in essentially the same manner as in [31J. Since u £ = v £ e 1 ^ and A £ = B £ + V^/io, we then 
immediately get the limiting dynamics for (u £ , A £ ) as well. 

Theorem 2 (proved later as Theorem 13.21) . Let j ex = h ex = 1 (parameter regime 1 ). Suppose that 
the initial data are well-prepared at order Co as defined by fll.32p . Then on any fixed time interval 
[0,T] the following hold. 

1. The vortex locations do not move in time, i.e. a»(t) = ctj(O) = for t G [0, T]. 

2. u £ — ^ weakly in if i 1 oc (fi\{a i } x [0, T]), where 

n / \ di 



M* = 

8=1 



x-_Oi\ e ^ +i/ _ ei0a+ ^ +i/) (136) 



f is defined by (11.281) . and ip* is a single-valued function on Q x [0, T] satisfying 

d t ip* — Aip* + ip* = V>*(0) in Q 
Vip* ■ v = — V6 a ■ v on dQ. 



;i.37) 



3. A £ — > A* in L 2 (Q x [0,T]) and h £ — ^ h* = curlA^, weakly in L 2 (Q x [0, T]]). The function h 
satisfies 

d t K - A/i* + /i* = 27r Yh=i d ^ai in Q ^ ^ 

h* = H ex on dQ. 
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Remark 1.1. If J = I = 0, then f — 0, and we recover a result from I31\/ . However, in the second 
item above, the equation satisfied by ip* is different from the equation satisfied by ip* in ISDj . The 
source of this disparity is the difference in choice of gauge. We work in a gauge where $ = /, 
whereas 131] utilizes the Lorentz gauge, $ + div A = 0. Formally changing (w*, A*, /) to the Lorentz 
gauge shows that our result is consistent with that of 131] . 

When 1 <C k ex <C | log £: | 1,/9 (regimes 2,3, and 4), the time scale X £ is much smaller than |loge|, 
and the vortices begin to move sooner. However, the limitation on the size of k ex still means that 
A £ — > oo so that it takes infinitely long for the vortices to begin moving in the limit e — » 0. In 
this case it is possible to extend the previous theorem to show that the vortices do not move in the 
limit, but unfortunately, the proofs of the second and third items of the theorem break down when 
1 <C k ex , so we can derive no information on the existence or structure of limits of v E , B £ , or curl5 e 
in the original time scale. 

The main result of the paper considers the solutions accelerated in time at scale X £ , i.e. we make 
the substitutions 

v £ (x,t) i — y v £ (x,X £ t) and B £ (x,t) i— > B £ (x,X £ t). (1.39) 

In this scaling the vortices move along well-defined, continuous trajectories. To show that the vortex 
trajectories are differentiable and to derive the limiting law governing their dynamics, we pass to 
the limit in a localized version of the evolution equation for the modified energy. We show that, as 
predicted, the applied boundary currents J ex and I ex exert Lorentz forces on the vortices in addition 
to the forcing term from the magnetic renormalized energy that was identified in [3Tj [26] . The exact 
form of the limiting law depends on the parameter regime. 

Theorem 3 (proved later in Lemma 14.11 and Theorem 14. 15j) . 

Suppose that the initial data are well-prepared at order Cq and that the solutions have been accelerated 
in time at scale X £ according to (11.391) . Suppose further than the initial vortex locations are separated 
from each other and the boundary by a distance at least a > 0. Then for < a* < a there exists 
a time T* = T^a*) G (0,T ] and n differentiable functions a { : [C^T*] — > Q satisfying the following. 

1. For each time t e [0, T*] there is a degree di vortex located at ai(t), i.e. the n initial vortices 
move along the trajectories ai. 

2. The vortices are separated from each other and the boundary by a distance at least a* for all 
time t G [0, T*] . In other words, the time T* is chosen to precede the first time at which a 
collision occurs or a vortex meets the boundary. 

3. If h ex = j ex = 1, then X £ = |loge| and the trajectories satisfy the dynamical law 

a t (t) = —VatWMt)) ~ 2^(V/ioK(t)) - V ± / (a,(t)) + V^/iKW)) (1-40) 

7T 

4- < h ex j ex <^ I log £: | 1/<9 ? then X £ = [log er| /j ex and the trajectories satisfy the dynamical 
law 

(H(t) = -2c/ i V ± A(a i (t)). (1.41) 

5- If0< j ex h ex <ti |log£:| 1//9 ; then X £ = |loge| /h ex and the trajectories satisfy the dynamical 
law 

Oi(t) = -2diVh (ai(t)) + 2d i V ± f (a l (t)). (1.42) 
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6. If I < h ex x j ex < |loge| 1/9 ; jex/Kx ->■ a, a™d h ex /k ex -)> /3 ; i/ien A £ = |loge| /fc e:c and i/ie 
trajectories satisfy the dynamical law 

d t (t) = -2d l aV ± f 1 (a l (t)) - 2d t f3 (V/i (a,(t)) - V i / Q (a i (t))) . (1.43) 

Remark 1.2. 

1. 5mce / zs completely determined by the applied boundary currents J ex and I ex , we see that 
the novel forcing terms in the dynamics, 2djV fk(a>i(t)) , k = 0, 1 are really forces induced 
by the applied currents. In this way we identify the new terms with the predicted Lorentz 
forces. The Lorentz force induced by the applied normal current, I ex , is always accompanied 
by a corresponding magnetic force 2diVh . On the other hand, the Lorentz force from J ex 
comes unaccompanied by a magnetic force, which indicates a fundamental difference between 
the currents J ex and I ex . 

2. In the case 1 <C k ex the magnetic renormalized energy disappears in the limiting dynamical 
law, and the motion of any given vortex is independent of the motion of the others. This is 
consistent with what was mentioned above about the form of the limiting dynamics for the 
gradient flow when h ex = 0(1) as compared to when h ex = (3 |loge|. 

3. In the third parameter regime, the current J ex dominates the dynamics and the limiting law is 
a perpendicular gradient flow. This implies that the dynamics of the limiting vortices possess 
n conserved quantities: Hi := /i(a.j) for i — 1, . . . , n. This flow keeps the vortices confined to 
the level sets of the function fx. Since we do not require f m J • v — (if J -v models a surface 
charge this should not hold), it is possible to choose J so that the level sets of fi form closed 
curves in Q, which gives rise to periodic motion of the vortices if no collisions occur. 

4- Working on the time interval [0, T*] prevents the vortices from colliding with each other or the 
boundary, and by letting <r* tend to 0, we can derive the dynamics up to the first collision time. 
Unfortunately, our techniques break down at a collision, and we can say nothing about what 
happens after. If a collision between two vortices occurs, it may be possible to use a deeper, 
more refined analysis as in J3, 0, [gpj/ to understand the dynamics afterward. On the other 
hand, if a vortex collides with the boundary, nothing in our analysis excludes the possibility of 
that vortex disappearing and another one nucleating somewhere else on dQ. 

5. Setting j ex = 1, J = I = 0, and H = 1, the second item of this theorem partially bridges the 
gap between what is known about the limiting dynamical law for the gradient flow dynamics 
when h ex = 0(1) l31f and when h ex = (3 |loge| l2bj . 

1.5 Plan of the paper 

The paper is organized as follows. In Section [2] we analyze the dynamics of the modified energy, 
proving it does not increase too quickly. In Section [3] we derive the limiting dynamics in the original 
time scale. In Section 0] we study the accelerated solutions and derive the limiting dynamics for the 
vortices, identifying the force induced by the applied boundary currents J ex and I ex . 

At the end of the paper we present two appendices. Appendix |A] contains a few remarks on the 
well-posedness and regularity of the equations (jl.6p with the new boundary condition V^w • u = 
iuJ ex ■ v as well as some a priori estimates needed in a couple places in the paper. We will also 
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need several results from the static analysis of the Ginzburg-Landau energy. These are collected in 
Appendix [Bl 
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2 Energy evolution 

2.1 Evolution of the modified energy 

The obvious starting point for an analysis of the behavior of a sequence of solutions (u E , A E , $ e ) to 
( II. 6p as e —7- is an examination of how the energy F £ (u, A) evolves in time. To understand this, 
we consider the Ginzburg-Landau free energy density 



A) := \ (\V A u\ 2 + i(l - \u\ 2 f + |curl A\ 2 ^j 



(2.1) 



Then d t F e (u, A) = J n d t g®(u, A), and we are led to consider the time derivative of the energy density. 
Lemma 2.1. For any triple (w,B,£) (not necessarily solutions), we have that 
d t g° £ (w, B) = div^w, V B w) + curl((curl B)(d t B + Vf)) 

- (dzw,A B w + ^(l- M 2 )) - (fltS + VO-(V ± curl S + (z«;,Vbw)). 

Proof. The result follows from a simple calculation and an application of the following lemma, which 
records the commutator relations for the covariant derivatives. □ 

Lemma 2.2. For any triple (w,B,£) (not necessarily solutions), the following commutation rela- 
tions hold for the covariant derivatives, df := dj — iBj, d% := d t + i£: 

d£d?w - dfd^w = iw curl B 
WbO^w -d(.V B w = iw(d t B + V£)- 

So, to understand d t F e (u,A), we apply Lemma EH] to (u, A, $) and integrate over Q. The 
integrals of the div and curl terms yield a non- vanishing boundary integral because of the conditions 
Vau ' v — iuJex ' v an d h = H ex . In particular, we find that for solutions to (11.61) 

d t F £ (u, A)+ I \dzul 2 + \E\ 2 = [ (iu, d$u)J ex ■ v - H ex E ■ t. (2.4) 
Jn Jon 

This boundary integral is inconvenient to work with, so we pursue an alternate strategy for studying 
the energy evolution. The idea is to modify u and A in a manner that turns the inhomogeneous 
boundary condition V^n ■ v = iuJ ex ■ v and h = H ex into homogeneous ones. This is accomplished 
by introducing the modified order parameter v = ue~ 1 ^ and the modified vector potential B = 
A — h ex 'V ± ho, where / and ho are the functions defined by ( 11.281) and ( 11. 24ft respectively. We 
emphasize that we are not making a gauge change since we do not make either of the changes 
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A i — y A -f - V/ or $ i— > $ — d t f . The choice of / and ho then imply that \7 B v • v — and curli? = 0. 
Throughout the rest of the paper we will write 

h' := curli? = curl(A — hexV^ho) = h — h ex ho. (2-5) 

The trade-off for switching to homogeneous boundary conditions is that the equations pick up 
forcing terms. We record the equation satisfied by (v, B, $) now. 

Lemma 2.3. Suppose the triple (u, A, $) solve (11. 6p . Let v = ue~ l * , B = A — h ex V ± ho- Then the 
triple (v, B, $) solve 



d t v + i$v = A B v + — (1 - \v\ 2 ) + 2iV B v ■ Z £ + ivf -v \Z £ \ 2 (2.6) 

and 



e 2 



8 t B + V$ = V ± ti + (iv, V B v) + \v\ 2 Vf - (M 2 - l)h ex V ± h (2.7) 
along with the homogeneous boundary conditions 



(\7 bV . u = ondttxR 
W = ondttxR 



(2-8) 



Here we have written h! = curll? and Z £ is defined by (1 1 . 3 j) . 

Proof. A direct calculation gives the PDEs. The boundary conditions follow via (ll.29p . □ 

In considering the evolution of the energy (in particular, when working with v, B) it is most 
convenient to work in a gauge for which the electric potential $ is set equal to the function / = 

jexfl — h ex fo. 

Lemma 2.4. We can fix a gauge so that $ = / and B(0) satisfies the Coulomb gauge, i.e div 5(0) = 
in tt and B(0) ■ v = on dtt. 

Proof. We make a gauge change v i-> ve 1 ^, Bi->B + V£, $ i— y := $ — dt£, where £ is given by 

£(x,t)= [ $(x,s)ds-tf(x)+ri(x) (2.9) 



and rj is the solution to 

'-An = divS(O) in tt , 

1 v ; 2.10 

Vn • v = -B(0) ■ v on dtt. 

Then \& = / and (5 + V£)(0) satisfies the Coulomb gauge as desired. □ 

Remark 2.5. In a// 0/ what follows we work exclusively in this gauge, which we call the $ = / 
gauge. We will also cease to refer to solution triples (v, B, $) and instead refer to just (v, B) since 
$ = /. 

Now we record the energy evolution equation for (v, B) in the $ = / gauge. 
Lemma 2.6. Let (v,B) solve f l23|) - fl2T8|) in the $ = / #a?j#e. T/ien 

d t g° E (v,B) = div(d t v,V B v) + cnr\(h'd t A) - \d t v\ 2 - \d t B\ 2 

+ {\v\ 2 - l)Z £ ■ d t B + 2(d t v,iV B v) ■ Z £ - (v, 8 t v) \Z £ \ 2 (2.11) 
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Proof. Plug the triple (v,B,0) into Lemma [2.11 and then plug in the equations of (I2.6p - fl2.8l) with 
$ = /. □ 

Since / and ho do not depend on time, the last term on the right hand side of the last equation 
is a time derivative. This leads us to define the modified energy density 

~g° £ (v,B):=g°(v,B) + ^\Z £ \ 2 (2.12) 

so that (v, B) satisfy 

d t [ g° £ (v,B) + [ \d t v\ 2 + \d t B\ 2 = [ (2(d t v,iV B v) + (\v\ 2 -l)d t B)-Z £ . (2.13) 
Jn Jn Jn 

The utility of this equation is that it allows us to identify the two distinct terms that contribute to 
the change of the modified energy. We find a standard dissipative term, 

\d t v\ 2 + \d t B\ 2 , (2.14) 

n 

which acts to decrease the modified energy. We also find the interaction term, 

/ (2(d t v,iV B v) + (\v\ 2 -l)d t B) Z £ , (2.15) 
Jn 

which mediates the interaction between (v, B) and the applied boundary currents J ex , I ex and mag- 
netic field H ex . We will eventually see that there is an exceptionally nice structure to the interaction 
term, but for now we ignore this structure and present a crude preliminary estimate of the growth 

o£j a s° e (v,B). 

Proposition 2.7. Let (v,B) solve f l2.6p - fl2.8p in the $ = / gauge. Then for any t>0, 

[ g° e (v,B)(t) + l [ [ \d t v\ 2 + \d t B\ 2 < exp(Ct) [ g° £ (v,B)(0), (2.16) 
Jn 1 Jo Jn Jn 

where C = 4 H^H^oo • 

Proof. Cauchy's inequality implies that 

/ 2{d t v,iV B v)-Z e < \ [ \d t v\ 2 + 4\\Z £ \\ 2 LOC f\\V B v\ 2 (2.17) 
Jn ^ Jn Jn * 



and that 



1 f 9 - o „ _ ,,9 f (1 — If 



(M z - l)d t B -Z £ <- I \d t B\ z + 2e z ||Z £ ||^ / v ^ )2 . (2.18) 
n 1 Jn Jn 45 

Plugging (I2.17P and (I2.18P into Lemma 12.131 then yields the differential inequality 

d t [ g° £ (v,B) + l [ \d t v\ 2 +\d t B\ 2 <A\\Z £ \\ 2 LOO [ g° £ (v,B). (2.19) 
Jn * Jn Jn 



An application of Gronwall then proves (12.161) . 



□ 
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2.2 Precursors for estimating B in the $ = / gauge 

In order to implement a more refined analysis of the interaction term (12. 15ft . we will require the 
ability to control the term J n \B\ 2 . Since curl B is gauge invariant and div I? can be controlled in 
some gauge choices, this is most naturally accomplished by using a Poincare inequality of the form 

B\ 2 <C I \divB\ 2 + \cm\B\ 2 . (2.20) 
Jn 

Such an inequality fails in general (e.g. B = Vh for h harmonic), but is available, for instance, in 
the space H^(Q; M 2 ) = {Be H 1 ^; M 2 ) \ B ■ v = on d£l}, which the magnetic potential B belongs 
to in both the Lorentz ($ + div A — 0) and Coulomb (div A = 0) gauges. Unfortunately, in the 
$ = / gauge, it no longer holds that B ■ v = 0, so we must resort to a version of Poincare that also 
involves a boundary term. In this section we will record such a Poincare inequality, then prove two 
lemmas needed to conveniently use this inequality in the $ = / gauge. 

We begin with the aforementioned version of the Poincare inequality, as well as another useful 
estimate. 

Proposition 2.8. Let 

X := {B e L 2 (tt;R 2 ) \ cui\B,divB e L 2 (tt), and B ■ v e L 1 \dVL)} . (2.21) 
Then there exists a constant C > so that for any B G X it holds that 

\B\ 2 < C ( I \cur\ B\ 2 + \div B\ 2 + I \B ■ u\ 2 ) . (2.22) 
'n \Jn Jan J 

Proof. For any B G X we solve the elliptic problem 



A0 = div B in 
V0 • v = B ■ v on dVt. 



(2.23) 



Then standard elliptic estimates give 

||0Li ( n) < C (j|div5|| i2(n) + \\B ■ Hlir-V2 ( an)) < C (j|div5|| i2(n) + \\B ■ u\\ L2(dQ ^ . (2.24) 

We may then define the vector field B = B — V<f), which satisfies div B = in Q and B ■ v = 
on dQ. The Poincare inequality mentioned at the beginning of the section is applicable to B and 
yields 

B <C curlS =(7||curlS|| ra , ov (2.25) 

L 2 (n) L 2 {n) n v l) 

Combining these two estimates yields the desired inequality. □ 

Remark 2.9. The Poincare inequality recorded in this proposition is not optimal, but is well-suited 
for our analysis. The above proof shows that it could be replaced with 

\\B\\ L2(n) < C (\\divB\\ L2(n) + ||curlB|| £a(n) + \\B ■ u\\ H . 1/2(dQ) ^ . (2.26) 

In order to apply this version of the Poincare inequality, we will need estimates for div B and 
B ■ v in the $ = / gauge. These are presented in the following lemma. 
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Lemma 2.10. Let (v,B) solve f l23|) - f[2~8j) in the $ = / #an#e. T/ien 

||S ■ v{t)\\ L , m < J \\{\v\ 2 - l) Jex Vf\ ■ v\\ L2{m) (2.27) 

and t t 

\\div B(t)\\ L2(Q) < ^ \\d t v\\ LH n) + I IKM 2 -!)/!!^) (2-28) 

Proof. The second equation in (I2.6l) - fl2.8p reads 

d t B = V ± ti + (iv, V B v) + (\v\ 2 - l)Z £ . (2.29) 

Taking the dot product of this equation with the boundary normal v and applying the boundary 
conditions h! = and V B v ■ v = on dQ yields 

d t B-v = V ± /i / • z/ + (iu, V B v • v) + (|u| 2 - 1)(V/ - h ex V ± h ) ■ v 

= (M 2 - l)(j ex J - h ex I + hexl) ■ V = (M 2 - l)j'exV/i ■ V. (2.30) 

Then, since B ■ z/(0) = 0, we get 

B ■ v(t) = [ (\v(s)\ 2 - l)j ex Vh ■ v, (2.31) 
Jo 

from which (12.27!) follows. 

Taking the divergence of the second Ginzburg-Landau equation, we find that 

d t div B + Af = div((iv, V B v) + \v\ 2 Z e ). (2.32) 

On the other hand, since Af = f, taking (iv, •) with the first equation in (I2.6p -( l2~8l) yields the 
equality 

(iv, d t v) + \v\ 2 f = div((iv, V B v) + \v\ 2 Z £ ) (2.33) 

so that 

Since div 5(0) = 0, we get 

div B{t) = I (iv,d t v) + I (\v\ 2 -l)f, (2.35) 



d t divB = (iv,d t v) + (\v\ 2 - l)f. (2.34) 



which yields (jZHP . □ 

The next result provides control of the boundary term (|t>| 2 — l)|an, which is necessary for the 
estimate (I2.27P to be useful. 

Lemma 2.11. Suppose v : Q — > C (not necessarily a solution) satisfies \\v \\ LOO < K and 

:|VH| 2 + (1 ~} V } )2 <fl~|loge|. (2.36) 



As* 



Then 



[ |1 - M 2 f < Cv^|log£| (2.37) 
Jon 

for a constant C > 0, depending on fl and K . 
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Proof. Using trace and interpolation theory, we have 



111 I l 2 ll 2 ^ r>\U I |2|| 2 s n\U I l 2 ll 3/2 111 I l 2 ll 1/2 /t ooN 

II 1 -M L 2(sn) <C||l-|^l || H 3/4 (n) <G\\-\v\ || ffl(n) ||l-H \\ LHQ y (2.38) 
Since \v\ < K, we may bound |V(1 — \v\ 2 )\ = |2 \v\ V \v\\ < 2K |V \v\\, so that (12.361) implies 

II 1 " \ v \ 2 \\m { n) ^ (^ 2 ^|VM| 2 + (1-M 2 ) 2 ) ' < V^llogel. (2.39) 
On the other hand, (12.361) also implies that 

||1 - M 2 || L2(Q) < V^Ke 2 \\oge\ (2.40) 
Then (12.371) follows by combining these three bounds. 

□ 



2.3 Control of the modified energy 

Fortunately, the crude Cauchy-Gronwall combination used in the proof of Proposition 12.71 is not 
the only method available for estimating the change of the modified energy. The second method 
relies on a more careful analysis of the interaction term ( 12.151) . We begin this line of reasoning by 
examining the space-time 1-form defined for any triple (w,B,^) by 

u := ((iw, dfiv) - g)dt + ((iw, d B w) + B), (2.41) 

where we have identified the vector B with the spatial 1-form B = B\dx\ + B 2 dx2 and written the 
1-form dsw := d sp w + iwB for d sp the spatial exterior derivative in Q. This 1-form is a modification 
of the usual space-time current, (iw,d^w)dt + (iw,dBiv). Taking the exterior derivative of u, we 
find the space-time Jacobian measure 

J = duo = Vidxi A dt + V 2 dx 2 A dt + /xofei A dx 2 , (2.42) 

where \i = curl((wu, Vbw) + B) is the usual spatial vorticity measure and the vector 

V := (V u V 2 ) = V(iw, d^w) - d t (iw, V B w) - (d t B + VO (2.43) 

is the "velocity component" of the space-time Jacobian measure. We will sometimes write fi = 
fjb(w,B) and V = V(w,B,£) to emphasize the dependence of fi and V on (w,B,£). Note, though, 
that both /j, and V are gauge invariant. 

Since d o d — 0, it holds that dj = 0, which implies the equation 

dtfi + cwlV = 0. (2.44) 

This relation allows us to identify the vector V with the velocity of the vortices. To see this, consider 
the simple example of /i(t) = 27r5 7 ( t ) for some smooth curve 7 : (io,£i) — > fi. A straightforward 
calculation then shows that V(t) = 2ii'y J -(t)5 1 (t). Note that even though V actually encodes the 
perpendicular to 7, we still refer to V as the "velocity." 

Now we relate the velocity associated to a triple (w,B,£) to a quantity that appears in the 
evolution equation. 
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2,n„.^ ( 2 - 46 ) 



Lemma 2.12. For any triple (w,B,£) it holds that 

V(w,B,£) + d t B + V£ = 2{d i w,iV B w) + \w\ 2 (d t B + Vf)- (2-45) 

Proof. Using Lemma I2.2[ we calculate 

W(iw,d^w) = {d^w.iV B w) + (w, Vb^w) 

= iV B w) + (iiu, 8^7 B w + iw(d t B + V£)) 
= (<%u;, iV B w) + d t {iw, V B w) - (id^w, V B w) + + V£) 

= 2(<%w, iV B w) + V B w) + \w\ 2 (d t B + Vf). 

Hence 

V(tu , 5, + (fi^B + VO = V(iw, dfiv) - dt(iw, V B w) 

= 2(d^w, iV B w) + \w\ 2 (d t B + VO- (2.47) 

□ 

We now use this lemma to rewrite the interaction term in the modified energy evolution equation. 
Lemma 2.13. Let (v,B) solve flZJ>])-(EE} in the $ = / gauge. Then 

d t ~g° £ (v, B) = div(d t v, W B v) + cur\(h'd t B) - \d t v\ 2 - \d t B\ 2 + V(v, B, 0) ■ Z £ (2.48) 

and 

d t [ g° £ (v,B)+ [ \d t v\ 2 + \d t B\ 2 = [ V(v,B,0).Z e . (2.49) 
Jn Jn Jn 

Proof. Use the triple (v, B, 0) in Lemma [2.121 to find 

(M 2 - l)d t B + 2(d t v, iV B v) = V(v,B,0). (2.50) 

Plugging this into Lemma [2.61 yields ( 12.4=81) . and ( 12.4=91) follows by integrating over Q and using the 
boundary conditions. 

□ 

The purpose of rewriting the interaction term as V(v, B, 0) ■ Z £ is that there are estimates of 
the velocity available from |25j. These estimates are actually for the velocity vector V(v,0,0), but 
they will be sufficient for analyzing V(v, B, 0). 

Proposition 2.14 (Product Estimate, Theorem 3 of [25]). Let u £ : Vt x [0,T] — > C satisfy 

r (i _ \ y iV 

/ |Vm £ | 2 + ^ 1 e| ; < C|logg| for all t G [0, T], (2.51) 

Jn 2e 



and 



Then the following hold. 



[ I |<W<C|loge|. (2-52) 
Jo Jn 
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1. There exist ft G L°°([0, T];M(Q)) and V = {V X ,V 2 ) with Vi G L 2 ([0, T];M(Q)) such that 
d t [i + cvx\V = and fx(u £ ,0) //, V(u e ,0,0) -» 7 as e zn x [0,T]))* /or any 

a G (0, 1). i/ere we have written Ai(Q) = (C°(Q))* for the space of bounded Radon measures. 



<liminf- ,( / / |Vn e -X| 2 / / \d t u £ \ 2 ) . (2.53) 



ti Jn Jti Jn 



2. For any [ii,^] Q IP>^1 anc ^ an ?/ e C°(f2 x [t 1 ,t 2 ];K 2 ), we nave i/ie bound 
i /■«. 

I V ■ X < I... 

e-s>0 | log £ | 

5. T/ie mapping t (//(£),£) is in if 1 ([0, T]) /or any £ G C£(0), and m particular 

Mh),0 - (fi(h),0\ < CVh^TilimM -=L= ( f 2 I \dtUsA 1 (2.54) 

^ V I log ^ I \Jh Jn J 

for any [ti, t 2 ] Q [0,T]. 



Remark 2.15. In the first item of Proposition 2.14 ^ e convergence we record is slightly stronger 



than what is stated in [25y . This is valid because the proof in [25y relies on the method used in \2_ 
where the stronger convergence result is actually proved. 

With this estimate of the velocity in hand, we can now show more refined control of the growth 
of the modified energy. The proof relies critically on the initial data satisfying a well-preparedness 
assumption given by (11.320 . Recall the dominant field strength, k ex , is given by (11.220 . 

Theorem 2.16. Let (v e ,B £ ) solve (12.6p - fl2.8p in the $ = / gauge. Suppose that the initial data 
(n e (0), B £ (0)) are well-prepared at order C , as defined by (ll.32p . Define the time scale \ £ := 
| log £r | /k ex . Then there exists a constant T > such that, as e — > 0, 

/ g° £ (v £ ,B £ )(t) < [ <7e(f e , -B e )(0) + 2C k ex for all t G [0, T A e ] (2.55) 
Jn Jn 



and 

,T A e 



o Jn 



\d t v £ \ 2 + \d t B £ \ 2 < 2C k ex . (2.56) 



Proof. The proof is inspired by that of Lemma III. 1 in [26]. For any t > 0, consider the two 
conditions 

I g° £ (v £ ,B £ )(s) < [ g° £ (v £ ,B £ )(0) + 2C k ex for all s G [0,t] (2.57) 
Jn Jn 

and 

ft f 

\d t v £ \ 2 + \d t B £ \ 2 < 2C k ex . (2.58) 



o Jn 
Define 

a £ := sup{t > | conditions (123TD and ( 12381) hold}. (2.59) 

Proposition 12.71 guarantees the existence of a time t £ > (depending on e, j ex , h ex and C ) such 
that both conditions hold for t £ . Hence a £ > for each e. We will show that actually a £ > T Q X £ for 
some T > as e — > 0, thereby proving the theorem. Suppose now, by way of contradiction, that 

liminf-^ = 0. (2.60) 

£^0 \ F 



We may suppose, up to extraction of a subsequence, that a e /X e — > as e — > 0. 

Rescale in time at scale a £ by defining w £ (x,t) = v £ (x,a £ t) and C s (x,t) = B e (x,a e t). By the 
definition of a £ , the inequalities 

/ g° £ (w £ , C £ )(t) < [ ~g° £ (w £ , C e )(0) + 2C A:« e for all t G [0, 1] (2.61) 
Jn Jn 



and 



— / / \d t w £ \ 2 + \d t C £ \ 2 < 2C k ex . (2.62) 
OCe Jo Jn 

both hold, but at time t — 1 one of the inequalities must be an equality since w £ and C £ are smooth. 
Our goal is to show that, using the product estimate, neither inequality can fail at time t — 1, 
producing the desired contradiction. In order to be able to apply the product estimate, though, 
we must first show that its hypotheses are satisfied. For the rescaled pair (w £ ,C £ ), equation (I2.49P 
becomes, after integrating in time from to 1: 



g° e (w e ,C E ){l)- I g° £ (w e ,C £ )(0) + - I [ \d t w £ \ 2 + \d t C £ \ 2 = [ [ V(w £ , C £ , 0) ■ Z £ . (2.63) 
n Jn a e Jo Jn Jo Jn 

Note also that for any time- independent vector field X G C°(Q; 1R 2 ), 



/ / V(w £ ,C £ ,0)-X= I [ V(w £ ,0,0)-X- [ [ d-, 
Jo Jn Jo Jn Jo Jn 



\{{l-\w £ \ 2 )C £ -X 

V(w e ,0,0)-X- I (1-K(1)| 2 )C £ (1)-X + / (l-K(0)| 2 )C e (0)-X. (2.64) 



'o Jn 

Proposition 12.101 implies that inequalities 



and 



/ \C £ (t) ■ u\ 2 < a 2 j 2 x ft (1 - \w £ \ 2 ) 2 (2.65) 

Jan Jo Jan 

[ \dWC £ (t)\ 2 <2 f [ \d t w £ \ 2 + 2a 2 \\f\\ 2 Loo f [ (l-\w £ \ 2 ) 2 (2.66) 
Jn Jo Jn Jo Jn 

both hold for all t G [0,1]. Condition (12.611) guarantees that the hypotheses of Lemma 12.111 are 

satisfied by w £ with a uniform constant K for all time t G [0, 1], so 

/ (l-K| 2 ) 2 <Cv^|loge| (2.67) 
Jan 

for a constant C > independent of t. Hence the inequalities 

/ \C £ (t) -u\ 2 <Ca 2 £ f ex \\V Ml^V^lloge] (2.68) 
Jan 

and 

/ |divC £ (t)| 2 <2a £ k ex + Ca 2 \\f\\ 2 Loc e 2 \\oge\ (2.69) 
Jn 

both hold for all t G [0, 1]. We also clearly have that 

/ |curlC e (t)| 2 < C|loge| (2.70) 
Jn 
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for all t G [0, 1]. We plug f l2T68|) -f l2T70|) into Proposition EH and recall that / = j ex f l - h ex f to 
deduce the bound 

/ \C £ (t)\ 2 < C(\\oge\ + a £ k ex ) + Ca 2 ^\hge\ k 2 x (\\f \\ 2 LOO + \\h\\ 2 LX + ||V/i|£„). (2.71) 
Jn 



Note that 



oieKx = -^■k ex X e = o(l) |loge| . (2.72) 

When combined with the trivial bound \Vw £ \ < \Vc £ w £ \ + \C £ \ and the bounds (I2.7ip . (I2.6ip . this 
implies that 

' f , Vw e \ 2 + {1 ~} W / )2 <C\loge\ (2.73) 



2V~" ' 2e* 

for all t G [0, 1]. This and condition (12.621) then allow us to apply Proposition 12. 141 to w £ to deduce 
the convergence of V(w e , 0, 0) to some V satisfying the bound (I2.53p . Plugging (12. 73 p . (I2.62p . and 
dekex I log e | into (I2.53P then yields 



o Jn 



V -X 



\x\ 



< 2 i-M^yc |loge| a £ k ex = o(l) (2.74) 
| log e | 



for any X G C°(Q x [0, 1]), and hence that V = 0. The vorticity measures /i(w £ , 0) also converge to 
a limiting measure /i, and the relation dtfi + curl V = then implies that d%\i = 0. We thus have 
that /i(t) = n(0) for all t G [0, 1]. Moreover, the bounds fl2UT|) and f l2TT|) imply that //(io e , C7 e )(t) -)- 
/i(t) = /i(0) for each t G [0, 1] as well, i.e. the vortices do not move. 

Returning to (]2.64j) . we employ the convergence V(w £ , 0, 0) — > and the bounds ( 12.61 j) . (I2.7ip 
to deduce that ^ 

[ [ V(w £} C £ ,0)-Z £ = o(l)k ex (2.75) 
Jo Jn 

Plugging this into (I2.63p . we see that 

/ ~g° £ (w £ ,C £ )(l) - [ ~ 9 ° £ (w £ ,C £ )(0) < o(l)k ex (2.76) 
Jn Jn 

We may therefore invoke Lemma IB.3I and Proposition IB. II to bound the modified energy at time 
t = 1 from below. Indeed, we find that 



/ g E (w e ,C e )(l) >7m|loge|+n7+- / \Z £ \ 2 + W d {a) + o{l)k e 
Jn 2 J n 

> [ g° £ (w £ , C e )(0) - C k ex + o{l)k ex . 
Jn 



(2.77) 



Inequality (I2.76P implies that (I2.6ip could not have been an equality at t = 1. Hence (I2.62p must 
have been an equality, i.e. 

— / [ \d t w £ \ 2 + \d t C £ \ 2 = 2C k ex . (2.78) 
a £ Jo Jn 

Plugging this and (I2.77P back into (12.631) then yields the inequality 

- C k ex + o(l)k ex + 2C k ex < o(l)k eX} (2.79) 

which is a contradiction as e — > since Co > 0. We deduce that it cannot be the case that 

liminf a £ /X £ = 0, i.e. that there exists a constant T > so that a £ > T X £ as e — > 0. □ 
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3 Limiting dynamics in the original time scale 

The following lemma allows us to remove singularities of the form {a{\ x [0, T] for solutions to 
parabolic equations with certain integral bounds. It is used in the subsequent theorem for deriving 
the limiting equation for the phase excess. 

Lemma 3.1. Suppose u is a weak solution to the equation dtu — Am = f on the set B(0,r)\{0} x 
[0, T] for some T, r > and f G L 2 (£>(0,r) x [0,T]). Further suppose that 

sup / K,, ( )|^ + /7 |V.«f,fa«<oo. (3.1) 

0<t<TjB{0,r) JO JB(p,r) 

Then u may be extended to a solution in all of 5(0, r) x [0, T]. 

Proof. It is shown in pQ that under the above assumptions, u may be extended to a solution of the 
equation on all of B(0, r) x [0, T] if and only if the thermal capacity of the set {0} x [0, T] vanishes. 
That this set has vanishing thermal capacity is established in Corollary 1 of [32]. □ 

Define a to be the multi-valued function, harmonic on f2\{<2j}, so that 



i=l 



Note that while a is multi- valued, its gradient is well defined away from the points {a{\. We are 
now able to derive the limiting dynamics in the original time scale. The following theorem, which 
is modeled on results in (211 I3T] . gives the dynamics for the pair (v E , B £ ), and since u e = v £ e^ and 
A £ = B £ + V ho, we may trivially derive the dynamics for (u £ , A £ ) from the theorem. 

Theorem 3.2. Suppose j ex = h ex = 1 so that J ex = J and H ex = H do not depend on e. Let 
(v e ,B e ) solve ( 12.61) - (12.81) in the $ = / gauge on Q x [0, T] for some T > 0. Further suppose that 
the initial data (v £ (0), B £ (0)) are well-prepared at order Co > in the sense of (11.321) . Then the 
following hold. 

1. /i e (t) — > yUo f or a H t G [0, T], i.e. the vortices do not move. 

2. v £ — v* weakly in f/" i 1 oc (f2\{a i } x [0,T]), where 

n / \ di 

Vm = Tl(±Z*) e ^ = e <e.-w. ( 3.3) 

' \X — Q,4' ' 



i=l 



where 9 a is defined by ( 13. 2ft and ^* zs a single-valued function on Q x [0, T] satisfying 

d t i/j* - A^* + ^* = ^*(0) «n ^ ^ ^ 

VV>* • v = — V6 a • z/ on <9fi. 

5. fi £ ->■ in L 2 (fi x [0,T]) and /i' e ^ = curl 5, wea% m L 2 (Q x [0,T]). Tae function h[ 
satisfies 

h[ = ondVt [ ' ' 
with the PDE satisfied in V'iVt x [0,T]). 
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Proof. We divide the proof into several steps. 
Step 1. 

We begin by showing that the energy stays well behaved and that the vortices do not move in 
time. Since the initial data are well-prepared at order Co we may apply Theorem 12. 161 to conclude 
that for e sufficiently small (so that T |loge| > T) we have the bounds 

/ g° £ (v £ , B £ )(t) < [ g° £ {v £ , B £ )(0) + 2C for all t G [0, T] (3.6) 
Jn Jn 

and 

/ / \d t v £ \ 2 + \d t B £ \ 2 <2C . (3.7) 
Jo Jn 

These bounds imply, as in the proof of Theorem 12. 16} that the limiting velocity vanishes, V = 0, 
and hence that fi £ (t) — > /xq for all t G [0,T]. This proves item[TJ 

Step 2. 

We now use the concentration of the energy around the vortex points {cij} to derive upper 
bounds for the energy away from the vortex points. Fix 

< a < ^min{dist(a;,<9fi)} U {{a* - aj\ \ % ^ j} U {6/V2} 

and define Q a = Q\ U" =1 B(ai,a). The bounds proved in Step 1 and the well-preparedness of the 
initial data allow us to apply Lemma IB. 31 to find that 

~ / |Vb^| 2 (*) + ^(1 - M 2 ) 2 (t) + I ! |curl5 £ (t)| 2 

< 7mlog- + n(7 + C) + W d {a) + 3C . (3.8) 
o 

for every t G [0, T}. 

From (I3.8P and Lemma f2.10[ we see that curl B £ and div B £ are both bounded when measured in 
the L°°([0, T]; L 2 (f2)) norm. On the other hand, Lemma [2.101 and Lemma [2.111 show that B £ ■ v — > 
in L°° ([0 , T}; L 2 (dQ)) . These facts, when combined with Proposition 12.81 and the elliptic Hodge 
estimate 

\\B\\ H i (n) < C (j|div£|| L2(n) + ||curl5|| L2(n) + \\B ■ v\\ H -x/2 (aa ^ (3.9) 

show that 

sup \\B £ \\ Hl(n) < (3.10) 

0<t<T 

for C* a constant that does not depend on e. On the other hand, ( 13.101) implies that d t B £ is 
uniformly bounded in L 2 ([0, T]; L 2 (Q)). This allows us to apply a result from [30] to deduce that, 
up to extraction, B £ — > B* in L 2 ([0, T]; H 1 ' 2 ^)). In particular, we also have that B £ — > B* in 
L 2 (fix[0,T]). 
Since 

|Vv 6 | < \^bM + \v £ B\ (3.11) 

we also deduce that up to extraction v £ — ^ v* weakly in H 1 ^^ x [0, T]) and v £ (t) — ^ weakly 
in /f 1 (f) (T ) for each t G [0, T]. Clearly v* is unit valued. 

Step 3. 
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We now further manipulate the bounds of the last step to derive the structure of v*. Since 
u* G -f/' 1 (fi cr ), we may write v* = e lLp * for <p* a multi- valued function such that V<^* is well-defined 
and Vv?* G L 2 (Q a ). Passing to the limit in ( 13. 8ft shows that 

J/ |V^(t)- J B,(t)| 2 <7 m log- + n( 7 + C , ) + ^ (i (a) + 3C (3.12) 

for every t G [0,T]. Define Y* = Vy* — V6 a , where 9 a is defined by (13. 2p . Then a simple 
modification of standard arguments (cf. [2T| [51]) shows that 

1 f \V<p*(t) - B*(t)\ 2 = ~ f \Y*(t)-B*(t)\ 2 + irn\og--0 (7 (l). (3.13) 

2 Jn a 2 Jn a ° 



Hence 



|y*(t) - 5,(t)| 2 < n( 7 + C) + W d (a) + 3C + 0,(1). (3.14) 



Since -B*(t) is bounded in L 2 (Q), we may let a — > in this inequality to conclude that Y*(t) is 
well-defined and bounded in L 2 {VL). 

The convergence // e — )■ 27r cfj5 aj , along with the convergence of v e and _B £ imply that 

n 

curl V</>* = 2tt c/j5 ai (3.15) 

in the sense of distributions. From this and the definition of O a we see that curl Y*(t) = in the 
sense of distributions for each t G [0, T\. By the weak form of the Poincare lemma (cf. Lemma 3 of 
[HI), it then holds that = V^>* for some 0* G if 1 (fi). Adjusting by a constant if necessary, for 
every t G [0, T] we may write if)*(t) = <p*(t) — a in Q a . Moreover, d t i/;* = d t <p* in x [0, T]. 

Step 4. 

We now derive the equation satisfied by ip*. Take (iv e , ■) with (I2.6l) - (l2.8p and expand: 

(iv £ , d t v £ ) + \v £ \ 2 f = div ((iv £ , Vv £ ) - \v £ \ 2 B £ + \v £ \ 2 Z £ ) . (3.16) 

For each a > the left side converges in T>'{Vt a x [0, T]) to d t ip* + f, and the right side to Aip* — 
div 5* + Af. Then, by (I1.28P Af = f, so we have that d t ip* — Aip* = — div B*. On the other 
hand, we know that d t divB £ = (iv £ ,d t v £ ) + (|f e | 2 — 1)/, so we may pass to the limit to find that 
d t div B* = d t ip*. Since div_B*(0) = 0, we find that div A* = ip* — ^*(0). We deduce that 

d t i>*- A^* + ^* = M Q ) ( 3 - 17 ) 

in n\{di} x [0,T]. 

It remains to show that the possible singularities at {a{\ are actually removable. To this end, 
we define the function 77* = ip* — '0*, where 

I I t/ a2 

Then 77* solves the equation d t r]* — A77* = ip*(0) — ip* ~ dt4>* G L 2 (Q x [0, T]). By the Poincare 
inequality 

r \Mt) -Mt)\ 2 < c [ |v^*(0l 2 <c* (3.18) 

Jn 
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for a constant = C(n, Q, h ex , d, a, Co). Hence 

sup / M 2 + f [ \Vr]*\ 2 < C*(l + T). (3.19) 

0<i<T JO 

We may therefore apply Lemma 13.11 to conclude that the singularities at {a{\ are removable and 
77* solves <9t?7* — A?]* = ~ i } * ~ dt^* in all of Q x [0,T]. This then implies that ip* solves 

(I3.18P in all of Q x [0,T]. The boundary condition Vy2* • v — carries over from the condition 
Vv £ ■ v = by multiplying ( 13. 161) by a test function that does not vanish on dfl but that vanishes 
in a neighborhood of the vortex locations and passing to the limit. So, V^* • v = — VO a • v on dQ. 

Step 5. 

We now derive the equation for h'* = curl B*. The magnetic potential B £ satisfies the equation 

d t B £ - V L curl B £ + B £ = {iv £ , V Bs v £ ) + B £ + (\v £ \ 2 - 1)Z £ . (3.20) 

We may take the curl of (13.201) in the sense of distributions and pass to the limit, employing the 
convergence of fi £ , to find that 

n 

d t K - + K = 2vr (3.21) 

i=l 

in V(Q x [0,T]). That h+ = on dQ carries over from the boundary condition h' £ = curl5 e = 0. 

□ 

4 Limiting dynamics in the accelerated time scale 
4.1 Preliminaries and vortex motion 

In this section we derive the limiting dynamics in the accelerated time scale. We rescale in time 
at the scale X £ := |loge| /k ex with k ex = max{h ex , j ex } by making the substitutions ( 11.391) . In this 
scaling the equations in the $ = / gauge become 



\ £ l d t v e = A Be v £ + ff(l - \v £ \ 2 ) + 2iV Be v £ -Z £ -v £ \Z £ \ 
X^dtB, = V L h' e + (iv e , V Be v £ ) + {\v £ \ 2 - 1)Z £ 



(4.1) 



along with the usual boundary and initial conditions. In the accelerated scale, the evolution equation 
for the modified free energy density becomes 

d t g° e {v e , B £ ) = div(d t v £ , V Be v £ ) + cm\(ti £ d t B £ ) - 1 \d t v £ \ 2 - 1 \d t B £ \ 2 + V(v e , B £ , 0) • Z £ . (4.2) 

We will assume throughout this section that the initial data (v e (0) , B e (0)) are well-prepared at 
order Co in the sense of (11.321) . Writing 

n 

Ai(0) = J>(0)^ (0) (4.3) 

i=i 

for the t = limiting vortex measure, we further assume that the initial vortex locations satisfy 

min{|ai(0) - aj{0)\ \ {dist(o f (0), dtt)} > a (4.4) 
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for some o"o > 0. We continue to assume that h ex ,j ex fall into one of the four regimes (|1.2ip . In 
any of the four cases A £ — > oo as e — > 0. 

Rescaling in time at scale X £ , Theorem 12. 161 provides for the existence of a constant T > such 
that 

/ g° £ (v £ , B £ )(t) < [ g° £ (v £ , B e )(0) + 2C k ex for all t G [0, T ] (4.5) 
Jn Jn 

and ^ 

[ ° / \d t v e \ 2 + \d t B £ \ 2 < 2C k ex X £ = 2C |loge| . (4.6) 
Jo Jn 

Note that the bound k ex <C |log 6r| 1//9 and (14. 6 p imply that d t B £ /X £ and d t v £ /X £ vanish in L 2 {Vt x 
[0,T ]), and up to the extraction of a subsequence we may assume that d t B £ {t)/X £ and d t v £ {t)/X £ 
vanish in L 2 (Q) for almost every t G [0, T ]. Recall also that from Lemma IA.21 the bound 
II ^b e v £ \\ Loo < C/e holds for all time. 

As the first order of business we record a lemma that shows the convergence of the space-time 
Jacobian and makes sense of the vortex trajectories. The result also establishes for any < cr* < o 
the existence of a time T* = T*(a*) G (0, T ] so that the vortex trajectories stay a distance ex* away 
from each other and dQ for all t G [0, T*]. We assume that such a a* is fixed throughout the section, 
and we work exclusively in the domain Q x [0, T*]. 

Lemma 4.1. Fix < cr* < cr . Then there exists aT t = T*(cr*) u>z£/i T* G (0, T ] so that 

1. The space-time Jacobian (ji(v e , B e ), V(v e , B E , 0)) ->■ (//, V) m (C ' a (O x [0,T ]))* /or a// a G 
(0,1). 

2. There exist functions G -ff 1 ([0, T*]; O), % — 1, . . . , n ; so i/iai for all t G [0, T#] 

min{|oj(t) - Oj-(t)| | i ^ j} U {dist(ai(*), SO)} > a* (4.7) 

and 

n 

H{t) = 27rX)d i (0) < y ai(t ). (4.8) 

i=l 

Proof. The bounds (14. 5 p and (14.61) allow us to argue as in the proof of Theorem 12.161 to deduce the 
first item. For the second item, Proposition III. 2 of [26J proves the existence of the vortex paths 
cij G -ff 1 ([0, T*]; fi), and the embedding if 1 ■=->■ C* ' 1 / 2 allows us to find the T* = T*{o~*) so that the 
vortex paths stay separated. □ 

Remark 4.2. Since the degree of the i th vortex does not change for any time in [0,T*], we may 
consolidate notation and write only di in place o/dj(0). 

According to Lemma |4.1[ the functions Oj G C 0,1 ' 2 ([0, T*]; 0) for z = l,...,n. So, for any 
< a < cr^/A, we may apply Lemma [B.3l to find that 



/ 



~g° £ (v £ ,B £ )<k ex C(a,n,tt) (4.9) 



and 



\h' £ (t)\ 2 <k ex C(a,n,n) (4.10) 
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for all t G [0, T*]. The latter bound implies that up to extraction 

Kit) 



V n ex 

whereas the former implies that 

Kit) 



Kit) weakly in L 2 (£l), (4.11) 



V ™ex 

We may then integrate both sides of the equation 

d t B £ 



Kit) weakly in H\Q\ U B(ai(t), a)). (4.12) 



A. 



+ (1 - \v £ \ z )Z £ = ^K + fa*, V Be v £ ), (4.13) 



against V -1 ^ for rj G C^°(Q), divide by ^k ex , integrate by parts, and pass to the limit to find that 
K(t) satisfies the PDE 

-AK(t) + K(t) = { 2 n ^ d ' S ^ "J-" 1 (4.14) 
[0 if k ex > 1 

in the sense of distributions. By trace theory and the fact that h' £ = on dQ, we also have that 
h'z = on dfl. Note that when k ex — 1, h'* is smooth outside of any neighborhood of U" =1 aj([0, T*]), 
and that when k ex 3> 1, K = identically. 



4.2 Convergence of the modified energy density 

The term in the energy evolution equation that allows for the identification of the vortex locations 
at each time is the normalized energy density g £ (v £ , B £ )(t)/ |loge|, viewed as a measure on Q. The 
mass of this measure is clearly bounded, so at any particular time we may extract a subsequence 
that converges in the weak sense of measures. The technical obstruction is that this extracted 
subsequence depends on the choice of t G [0, T*], whereas we would like the subsequence to converge 
for all t G [0, T*]. Fortunately, the measures satisfy a certain semi-decreasing (in time) property 
that allows us to find such a subsequence by adapting results from [5]. We begin with a proof of 
this semi-decreasing result. 

Lemma 4.3. Let (f> G C 1 (fi). Then for any < t 1 < t 2 < T* it holds that 

2 g° £ (v £ ,B £ )(t 2 ) f 2 g° e (v £ ,B £ ) ih) , , r , , , , 

1 ■ < C(n,C ,f,h ,(f))Vt2 -t x + o(l,0), (4.15) 



|loge| Jn |loge| 

where C(n, C , f, h , (ft, Q) is a constant depending on n, C , || | V/i| + | V/o| + | v7iq||| loo , and \\<fi\\ c i , 
and where o(l, <fi) — > as e — > 0, with rate of convergence dependent on </> but independent ofti,t 2 . 

Proof. Lemma 12.131 rescaled to the accelerated time scale and then integrated in time from t\ to 
t 2 , shows that 



k2 g° £ (v £ ,B £ )(t 2 ) f -g^B^t, 



loge| J n I log e | A e |lo. 



I + 11 (4.16) 
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where t 

/ := T7~~ r / 2 / ■ (d t v £ , V Bt v e ) + 20V ± • h' £ d t B £ (4.17) 

and ( 

// := — / 2 / 2 V(t; £ , fl e , 0) • Z e . (4.18) 
I log £ | Jn 

We apply Cauchy-Schwarz and bounds (14.51) and (14. 6 p to / to derive the inequality 

,„ < Wc , ( ti m y /2 ^ r/iw + i^/ i/2 



ii 



l°g £ l / \|loge|y tl 



<C(n,C o ,0)v / t 2 ~^i'. (4.19) 
To handle JJ first note that, according to Lemma 14. 1[ for a fixed vector field X G C°(Q; M?) 

2 / <P 2 X ■ V(v e , B e , 0) = / 2 / 0V • X + o(l, 0), (4.20) 
ti Jti Jn 

where o(l, 0) — » as £ — >• 0, with rate of convergence dependent on cf> but independent of t\, t 2 . We 
then apply item 2 of Proposition 12.141 along with the bounds (14. 5 ft and ( 14. 6 ft to bound 



' <C{n,C Q ,^X).y/t^h. (4.21) 



</> 2 X- V 
ti Jn 



Applying this with X = V/i, V/o, V^/iq anc l employing the fact that fc e:r <C |loge| then shows that 



|//| <o(l)Vh^ti + o(l,(j)). (4.22) 
Then (fl~T5D follows by plugging the bounds (gjgj) and (fl~22D into (|QS]i . □ 

With this semi-decreasing property established, we can now show that up to the extraction 
of a single subsequence, the modified energy density converges to a sum of Dirac masses for all 
te[0,T„]. 



Proposition 4.4. There exists a subsequence so that 

9°e(v £ ,B £ )(t) 



I log £ I 



n 

_ 

1 



weakly-* in (C 1 ^))* /or a// t e [0,TJ 



Proof. We will first establish the convergence of the measures and then establish the structure of 
the limiting measure. In order to establish the convergence (up to extraction) for all time, we will 
use a variant of Helly's selection principle in conjunction with the semi- decreasing property proved 
in Lemma [4.31 We essentially follow the strategy presented in Section 5.4 of [5]. 

Let {4>k}kLi be a countable set of functions in C 1 (fi) so that the span of {c/) 2 ,}^ is dense. For 
each k,e define the function £ fcj£ : [0, T#] — > R by 



~g»(v £ ,B £ )(t) 



m*) = ,r ; w • (4.24) 

I log e | 
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By Lemma 14.31 the functions ^f.,e satisfy the following semi- decreasing property: for every 5 > 
there exist > and > so that for every t 2 G (0, T) and t\ G (t 2 — Tk, t 2 ) it holds that 

&M(f 2 ) < €fc >e (*i) + for all £ < e fc . (4.25) 

Then a semi- decreasing variant of Helly's selection theorem (cf. Lemma 5.4 in [5]) implies that 
there exists a set of functions : [0,T*] — >■ 1R such that up to extraction 

£fc, e 0) -» for all t G [0, TJ and for all fc G N (4.26) 

as e — )■ 0. From the density of the span of we deduce that there exists a family of measures 

u(t) such that 

A weakly- * in (C 1 ^))* for all t G [0,2;]. (4.27) 

We now derive the structure of the limiting measures f(t). According to Lemmas 14.11 and IB. 3 

for a < er*/4 

~9° e (v £ ,B £ )(t) 



n\UB(ai(t),a) 



\loge\ 



< o(l). (4.28) 



Since a can be taken to be arbitrarily small, this then implies that 

n 

'(t) = X)o i (t)tf ai(t) . (4.29) 



i=l 



We now calculate the value of otj(t). Fixing a < a*/4, for each z we may choose r)i G C 1 (fi) so that 
supp(r/j) = B(a>i(t), 2a), rji — 1 on 2?(cij(i), a), and < 77* < 1. Then applying Lemma IB~2l with this 
choice of a, we have that 



<t) = lim f > lim inf / «?(;-.■ *.)(*) > „ (4 30) 



'B(ai(t),(T) 

On the other hand, from the bound (14.51) it holds that 

> aAt) = lim > / r\i : < lim sup / : < Tin. 4.31 

Hence aj(£) = 7r for each i — 1, . . . , n. 

□ 

Remark 4.5. VKe assume in what follows that we are working with the extracted subsequence so 



that the convergence result of Proposition 4-4 holds 



4.3 Convergence of the stress-energy tensor 

The stress-energy tensor associated to a configuration (v £ ,B £ ) is the symmetric 2-tensor, T £ , with 
components 

(T £ ) tJ = (d*°v £ ,d?*v £ ) - 5 tJ f 1 - \V Be v £ \ 2 + ±(1 - \v £ \ 2 ) 2 - i(curl5 e ) 2 ) , (4.32) 
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where dfv := dfl — iBiV is the i th covariant partial derivative. The divergence of T £ is the vector 
divT £ with components 

(diyT e ) i ^d 1 T li + d 2 T 2i . (4.33) 

The divergence encodes the "force" acting on the vortices, and by passing to the limit in T £ for 
solutions (v e , B £ ) we will be able to derive one of the terms driving the limiting vortex dynamics. 
The method of passing to the limit in the stress-energy tensor has been used extensively in the 
study of Ginzburg-Landau dynamics PU I2TJ E2J EH ESI E] ■ 

To make sense of the limit of T £ , we follow a strategy similar to that of Chapter 13 of |27j, 
where they study the limit of T £ for (v e , B £ ) solutions to the elliptic Ginzburg-Landau equations 
on Q (critical points of the Ginzburg-Landau energy functional). Roughly, the idea is to use the 
vanishing of the right hand side of the equation 

V x /£ + (iv £ , V Be v e ) = ^ + (1 - \v e \ 2 )Z e 

to show that T £ has the same limit as a similar tensor defined in terms of the induced magnetic 
field h £ , and then to derive the structure of the limit of the latter tensor. In [27], they prove a 
convergence in finite-parts result for T £ on all of Q, which is stronger than what we shall prove here 
for Q x [0, T*]. Indeed, the convergence result we prove here holds only for cylinders U r x 1 2 ] for 
which we know certain strong bounds on the energy. Here we have written U r for a ball of radius r 
rather than B r to avoid confusion with the vector potential B e . 

Theorem 4.6. Suppose that in a ball of radius r > 0, U r C Q, it holds that 

[ g e (v e ,B e ){t) < Ck ex for all t E [t u t 2 ]. (4.34) 

JUr 

For the limiting induced magnetic field, h'^, define the tensor 

S(K) := -VK ® V/< + / 2x2 (\ \VK\ 2 + 1 \h'A . (4.35) 



Then for < S\ < r, 



^-T £ -> S(K) in L\U S1 x [UM])- (4.36) 



Remark 4.7. Recall that if k ex ^> \, then the limiting magnetic field vanishes, i.e. h'^ = 0. So, in 
this case, the theorem says that T e /k ex — > in L 1 (U Sl x [ti,^])- 

The proof of the theorem is based on the following three lemmas, all of which rely heavily on 
the energy bound (14.341) . We begin by showing that \v £ \ must be close to 1. The argument is a 
modification of one used in 



Lemma 4.8. Fix < si < s 2 < r. Then for almost every t E [t\,t 2 ] it holds that 

Kmtf* \\1 - \v&)\\\l-w S2 ) = ° ( 4 - 37 ) 



e^O 



Proof. Fix t E [ti,t 2 ] so that d t v £ (t)/X £ -> 0, d t B £ (t)/X £ -> in L 2 {9). We know that the set of such 
t has full measure. For the rest of the proof we will neglect to write the dependence of the functions 
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on time, but all are implicitly evaluated at the chosen time t. We assume that e < (r — s 2 ) /^L so 
that k exy /e) C C/ r for all x G C/ S2 . Suppose, by way of contradiction, that 



lim sup A* 7 



£->0 



lL°°(t/ S2 ) 



> 0. 



We may extract a subsequence (still denoted by e) so that 



Ay 8 ill - 



\L°°(Us 2 ) 



> a > for all e. 



(4.38) 



(4.39) 



Since each v £ is continuous, we may choose points x £ G U s , 2 with the property that |t> £ (x e )| < 

1 \-V8 

1 — ctA e 

Rewrite the equations (14. II) in elliptic form as 



A ft v e + ^(l-|v e |" 
V x /i' + (iv £ , V Be v £ ) 



= K £ 
X 



where 



d t v £ 
A, 



2zVr v e ■ Z £ + v £ \Z e \ and X £ :-- 



d t B £ 



+ (i - \v £ nz £ . 



(4.40) 



(4.41) 



Since the quantities of interest are gauge invariant, we are free to switch to the Coulomb gauge 
(again, only at the time t) via v £ w £ := v £ e %e - e , B £ \-t C £ := B £ + V£ £ with £ e chosen so that 



div C £ = in U r 
C r ■ v = on 



(4.42) 



Then the elliptic equations satisfied by (w e , C £ ) are 

'A CE w £ + ^(l-\w £ 



V^K + (iw s , V Ce w £) 



(4.43) 



where k £ , X £ are as above (still with v £ in their definition) and h' £ = curli? e = curlC e . In the 
Coulomb gauge, we have that (cf. Proposition 3.3 of [27] ) 



\\r lr 

which implies that \\C £ 
Note that 



<C\\h' e r Hl{Ur) <C\ l U%'_..B;) 

< Ck ex . 

kexy/e 



d t B £ 



:i - \v. 



|2\2 



(4.44) 



\L°°(Ur) 



/ ~g° £ (w £ ,C £ )> / _~g° £ (w £ ,C £ )= / Is 



9° £ (w £ ,C £ 



dB(x e ,s) 



ds 



\loge\ 



which implies that 
inf 

k ex e<s<kexVe \ JdB(x E ,s) 

so that there exists r £ G (k ex e, k exy /e) with the property that 

AC 

~g° £ (w £ ,C £ ) < — 

dB(x E ,r e ) A e 



~g £ (w £ ,c £ )< 



I log £T | A e 



(4.45) 



(4.46) 



(4.47) 
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We now recall the Pohozaev identity on B(x £ ,r £ ) for solutions to (14.4=3)) : 

/ ^(1 - M 2 ) 2 - (K) 2 = [ ({e*>K„ V Ce w E ) - h! E Xi) ■ (x - x e ) 

JB(x E ,r E ) Z£ JB(x E ,r E ) 



+ r £ I \ \V CE w e ■ r| 2 - \ \V Cs w £ • u\ 2 + -1(1 - \w £ \ 2 ) 2 - \{h' £ f. (4.48) 

JdB(x E ,r E ) A A 4£ A 

We can use the L 2 (U r ) bounds on K £ and X £ that come from (14.34)) to control the first term on the 
right side, and ( 14.47)) controls the second term: 

f 1 AC AC 

/ ^(1 " \w £ \ 2 ) 2 ~ {Kf < Ckl x r £ + < CA£VS+ ^. (4.49) 



Now, by Sobolev 



/ (^) 2 < Cr E \\h' £ \\ 2 LHUr) < Ck ex ^e \\h'£ mUr) < Ck^y/i. (4.50) 

JB(x E ,r E ) 



(4.51) 



'B(a;e,r e ) 

Hence 

/ ^(l-H\ 2 r<Ckt x V-e+f<y 

JB(x E ,r E ) A£ A e *e 

for e sufficiently small. 

On the other hand, we have that 

\Vw £ \ < \V Ce w £ \ + \C £ \ = \V Be v £ \ + \C £ \ <- + Ck ex < ^. (4.52) 

e e 

Let s £ = a(2Co)~ 1 £:A i r 1//8 , and note that s £ < r £ for e sufficiently small. Then the mean value 
theorem shows that 

-^jg < 1 - \w B (x)\ 2 for all x e B(x £ , s £ ), (4.53) 

and hence 

vra 4 r (i - \w. |2 ^ 2 



32C 2 A £ 1/2 Jb(x.,b.) 2 ^ 2 
Comparing (I4.51f) and ( 14.541) . we deduce that 

net 4 C 
32C 2 AP V 

which yields a contradiction as £ — > 0. 

□ 

The next lemma, which is a parabolic modification of the elliptic result in Proposition 13.4 of 
[27] , shows that we can essentially replace V Be v £ with (iv e , V Bs v £ ) in the definition of T £ . Note that 
this lemma is the only place we use the full strength of the upper bound k ex <C llogej 1 ^ 9 . 

Lemma 4.9. Define the tensor 

T' £ = (iv £ ,V Bs v £ ) ® (iv £ , V Bs v £ ) - J 2x2 Q \(iv £ , V Be v £ )\ 2 - \(K) 2 ^ . (4.56) 

Then 

k e x ll^e _ T £ \\L 1 (U ai x[tiM) ~^ ®' (4-57) 
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Proof. By Lemma 14.81 for almost every t G [ii, £ 2 ], as e — >• 0, we may write v £ (t) = p £ (t)e iVe ® for 
p £ , <p £ well-defined and single- valued in U S2 . We again neglect to write the dependence on t in what 
follows. Then since (iv £ , Vb s v £ ) = p^(V<p £ — B E ), 



P % - t> = p\ ( v Pe ® v Pe - / 2X2 ( \ |v Pe | 2 + ) + - 1) { - 2 



'ef 



(4.58) 



We may then write 
so that 



t £ - r e = (i - p 2 £ )T £ + p 2 £ t £ - r £ 



The first term on the right is easy to manage in view of (14.341) and Lemma 14.81 



I? 



C 



[ \l-p 2 \~g° £ (v £ ,B £ )^0. 



(4.59) 
(4.60) 

(4.61) 



To handle the second term, we take (v £ , •) with the first equation of (14. ip to find that p £ satisfies 
the PDE 



dtps 
A £ 



- Ap £ + p £ |V<p £ + Z £ - B £ \ 2 = %(1 - p 2 £ ). 



(4.62) 



Multiply (I4.62p by (1 — p £ ) and integrate over a ball U s for < s < s 2 to be chosen later. After 
integrating by parts and rearranging, we arrive at the equation 



/ | V p £ | 2 + ^(l-p 2 )(l-p £ )= / p £ (l-p £ )|V^ £ + Z £ -B £ 



dtps 



+ I -f^(l-p 6 )- 

<U S A e 



lt (i -*> (4 - 63) 



By a mean- value argument, we may choose s G [si, s 2 ] so that 



9 2 

|Vp £ | 2 < , 

IdUs S 2 — Si J Ur 

On U S2 , for e sufficiently small, it holds that 



9°e(Ve,B £ )< 



Ck e 



S2 - Si 



(4.64) 



4 " l + p £ 
Combining (I4.63[) ( 14.65[) . we find that 



i<-*- and ^< ° 



(4.65) 



1 f 1 ,„ |2 (1-P £ ) C 



4e 2 



k P _ x X £ Ju : 



\^B e v £ + iv £ Z £ 



C 



d t v £ 



A, 



1/2 



v/2 



TVS 



Ckp X \ Ck e 



S 2 - Si / A, 



o(l), (4.66) 



where the last equality follows from the fact that A £ = [log er| /k ex and k ex <C |loge 



1/9 
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Hence 

— / \T £ - T' s \ (t) ->■ for almost every t G [t u t 2 \. (4.67) 

hex JU S1 

Since for all t G [^1,^2], 

~r~ f \Te~ T' £ \ S e ) < (4-68) 

K ez JU S1 K ex JUr 

we conclude that (14. 57ft holds via an application of the dominated convergence theorem. 

□ 

The third lemma establishes the strong convergence of h' £ /y/ k ex and Vh' £ / \Jk ex via an argument 
like that of Claim 5.4 in EH. 



Lemma 4.10. It holds that 

hi 



— ->■ K and —j= -> V< m L 2 (t/ Sl x t 2 ])- 



Proof. Since V" 1 -/^ = — (iv e , Vb s v £ ) + d t B £ /\ e + (1 — \v e \ 2 )Z e , we know that 



(4.69) 



v ex JU, 



Wh'f < — 



v B v £ \ 2 + (i-\v £ ?y\z £ \ 2 + 



d t B £ 



X, 



(4.70) 



so for almost every t G [ti, t 2 ] it holds that 

±- [ \h'f+\Vh' £ \ 2 (t)<C. 

K ex Ju r 



(4.71) 



Then h' £ /y/k^ — ^ h* weakly in H l {U r ), h!J\fk^ x —> h'^ in L 2 (U r ), and by dominated convergence 
K/VKx -»■ ft,* in L 2 ([/ s x [t 1 ,t 2 ]) as well. Write v £ = p £ e %ipe in f/ S2 . Then 



1 d t B £ 1-p 



P 2 e A £ p £ 



2 ±Z e -V<p e + B e = ^V x h' e . 



1 



(4.72) 



Fix £ G C%°(U S2 ). Multiply the last equation by V _L (£(^ — ft*)) and integrate over U r to get 



U r P . 



^(IVft'j -vftt-vftl 



C/r Pe 



+ 



' ' ^Z £ )-Y^Uh>-h>)). (1.73) 



p2 



From this and the known convergence results, we deduce that for almost every t G [£1,^2 



and hence that 



^ f f |vft £ | 2 (t)^/ ZWK\ 2 (t) 

Kex Ju r Pe JU r 

J_ /" £| V ft4-Vft'J 2 (t)^0. 

K ex JUr 



(4.74) 



(4.75) 
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Now take £ so that < £ < 1 and £ = 1 on U Sl to see that 

— / \Vh' e - Vh'f (t) ->■ for a.e. t G [*i,t 2 ]. 



(4.76) 



We cannot apply the standard dominated convergence theorem directly since we lack a function 
that dominates the term J v \d t B £ /X £ \ 2 in (14. 70 p . However, 

2 



ex JU, 



iV/i'J < c 



c 



ex JU, 



d t B e 



X, 



and 



ti 



ex JUr 



d t B £ 



— >■ 



c, 



so by a variant of the dominated convergence theorem (cf. section 1.3 of [i~5 



|V^-V/i'J ->0. 



(4.77) 
(4.78) 

(4.79) 
□ 



We are now in a position to present the 
Proof of Theorem \4-6] We have that 



d t B £ 



A, 



+ a-kr)z e ->o 



in L 2 (U r x [ti, ^2])- When combined with Lemmas 14.91 and I4.10[ this proves that 

K ex \z, Z 

in L 1 (U S1 x [t 1? t 2 ]). The result follows by noting that for any function £ : IR 2 — > 



(4.80) 

(4.81) 

(4.82) 
□ 



We conclude this section with a lemma that links the tensor S to the renormalized energy in 
the case j ex = h ex = 1. The result is an adaptation of Application 3 in Chapter 13 of [27] . 

Lemma 4.11. Let h solve 

-Ah + h = 2n J27=i di?>ai in ^ 
h = on dVl 



(4.83) 



and define the tensor 



S(h) = -Vh ®Vh + J 2x2 ( \ \Vh\ 2 + I \h\ 2 



(4.84) 



(4.85) 



Then for r > sufficiently small so that a,j ^ B(a,i, r) for j ^ i, it holds that 

[ S(h)v = V ai W d (a) + o r (l), 

JdB(ai,r) 

where d = {d±, . . . , d n ), a = (a\, . . . , a n ), Wd{a) is the renormalized energy defined by ( 11.121) . and 
o r (l) means a quantity vanishing as r — > 0. 
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Proof. We begin by expanding S(h)u in the orthonormal basis (v,t) on dB(a>i,r). This yields 

S(h)u = - ((Vh ■ t) 2 - (Vh ■ u) 2 + h 2 )v+ (Vh ■ u)(Vh ■ r)r. (4.86) 



2 

Now we use the decomposition 



it 

h = J^djiM-^i) - log I' " a M ( 4 - 87 ) 



where Sq is defined by ( 11.131) . This allows us to write 



Vh(x) = J2 d A V^(x,a,) - X a \ ) =: -di X a \ + VHi(x). (4.88) 



Note that the function Hi E C l > 1,2 (B(a h r)). Then on dB(at , r), we have that 



Vh • v = -— + VH { ■ v and Vh ■ r = VH { ■ r, (4.89) 
r 



which implies that 



= U~$ + ^tA » + (™^)t + 0{ 1). (4.90) 



Integrating, we get 

I S(h)v=-f VHi + o r (l) =2irVH i (a i ) + o r (l). (4.91) 

JdB{ai,r) r JdB(a t) r) 

It is then straightforward to check that VHi = V ai Wd(a), and the conclusion follows. □ 
4.4 Dynamical law 

We have now established all of the preliminary convergence results necessary to derive the dynamical 
law. As the first order of business, we calculate the divergence of the stress-energy tensor, T £ . 

Lemma 4.12. Let T £ be the stress- energy tensor associated to (v £ ,B £ ), defined by (I4.32p . and let 
divT e be its divergence vector, as defined by ( 14.331) . Then 

divT £ = 1 ((d t v £ , V Bs v E ) - h' £ d t B^) + (v £ ,V Bs v £ ) \Z £ \ 2 -ti(v e ,B e )Z^ (4.92) 

Proof. A straightforward calculation shows that divT e = (K e , V Be v £ ) — h' e X £ , where 

K £ :=^- 2iV Be v £ -Z £ + v £ \Z £ \ 2 and X e := ^ + (1 - \v £ f)Z £ . (4.93) 

The result follows by plugging K £ ,X £ into the formula for divT e and noting that for any Zsl 2 , 

{-2iV Be v £ ■ Z, V Be v £ ) = 2{d?°v e , id^v £ )Z L (4.94) 

and that 

(\v £ \ 2 - l)h' £ + 2(d^v £ , id B °v £ ) = - cux\(iv £ , V Be v £ ) - h' £ = -fi(v £ , B £ ). (4.95) 

□ 
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We now record a local version of the energy evolution equation that contains the divergence of 
the stress-energy tensor. 

Lemma 4.13. Fix <\> G C£°(fi). Then 

d t f i 9 °j^ Be) + — J—, / (|^ e | 2 + \d t B £ \ 2 ) = [ <PV(v £ , B £ , 0) ■ Z £ 

J n |loge| A e |log£-|./ n XekexJa 

2 



/ V0- (divr e -vK| 2 i^L + M« 6> 5 e )^) . (4.96) 



Proof. Note that 

div(<9 t ^, V Be?)£ ) + curl(/i' e a t S £ ) = div((a^ £ , V S ^ E ) - ^5^). (4.97) 
Using this, we rewrite ( 14. 2 p using Lemma [4.12| multiply by <j)/ |loge|, and integrate by parts. □ 

The following result is the main ingredient in deriving the dynamical law for the vortices. It 
combines all of the previous convergence results in the case when we know that a vortex path is 
contained in a given cylinder and the energy is bounded in a wider cylinder. 

Proposition 4.14. Suppose dj(t) G B(xo,ri) for all t G [ti,^] and that 

! g° £ (v £ ,B £ )(t)<Ck ex (4.98) 

J B(xo,r2)\B(xo,ri) 

for all t G [ti, ^2] ■ Then for any <fi G C™(Q) with the property that 



supp(L> 2 0) C {x G Q I r\ < \x — Xq\ < r 2 } 

it holds that 

<<l>Mt 2 )) - tMh))) 

rt 2 r r*2 



rt2 

D 2 <p : S(K) - 2ixd 3 I Z ± {a j (t))-V4>{a j {t))dt, (4.99) 
Jti 



't\ J B(xo,r2)\B(xo,ri) J t\ 

•where Z : Q — > R 2 is given by 

Z = lim-A (4.100) 

6^0 k ex 

Proof. We integrate the result of Lemma [4. 131 from t\ to t 2 to find 

I + II = III + IV + V (4.101) 

for 



I : 



.g° e (v £ ,B e )(t 2 ) f J° £ (v £ ,B £ )(t 



1 j 



I log £ I J n |loge| 

//:=t-^— r r / <P {\d t v £ \ 2 + \d t B £ \ 2 ) , 
A £ I log e I J tl J Q 



12 



ti J ft ^e^ex 



III:= I I ^—V(v £ ,B £ ,0)-Z £ 

z^ex 
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hex jt\ Jn 
and 

r-ta 



TV := f 2 [ V0- (divT £ - V|w e | 2 — ' 



V:=- I' [ vi{v e ,B g )±Z±.V<l>. 

Jti JQ K ex 



We will pass to the limit e — > in each term. 
From Proposition I4.4[ we know that 

/^ 7r (0( Gj .( t2 ))_0( a .( tl ))). (4.102) 

The bound (14. 6p . together with A e — > oo implies that JJ — )■ 0. The convergence of V(u e , 5 e , 0) and 
the multiplication by 1/A e show that J/J — >■ 0. To handle IV we first note that 

i / V0-VK| 2 ^ = -^ / v0-v(K| 2 -i)|z £ | 2 



^- / (1 - K| 2 ) div (V0 |Z £ | 2 ) = o(l). (4.103) 



1k P 

By assumption, the support of D 2 <p is contained in the interior of B(xo,r 2 )\B(x ,ri), and hence 

2 / V0-divT £ = /" 2 /" D 2 <P:T £ (4.104) 

il ^ J B(x ,r 2 )\B(x ,ri) 

Using the energy bound (14.981) . Theorem 14.61 and a covering argument, we deduce that 

-> in L 1 (supp(D 2 0) x [t 1} t 2 ]). (4.105) 

"'ea; 

Hence ^ 

IV ^ f 2 [ D 2 <P:S(K). (4.106) 
Jt\ Jn 

For V we use the convergence /j,(v e ,B E ) — > 2nJ2diS ai (t) in conjunction with item 3 of Proposition 
EH to get 

V -> -2^^ / Z- L (a i (t)) • V0(a i (t))rft. (4.107) 
Jti 

□ 

We now derive the dynamical law by using this result with an appropriate test function. 
Theorem 4.15. The vortex trajectories Oj(t), i = 1, . . . ,n, are differentiable. 

1. If hex — J ex = 1; then the trajectories satisfy the dynamical law 

Oi(t) = --V a ,W d (a(t)) ~ 2d i (Vh (a i (t)) - V ± f (a l (t)) + V ± f 1 (a l (t))) (4.108) 

7T 

2. In the other three parameter regimes fll.21f) we define 

lim ^ = a E [0, oo) and lim = e [o, oo). (4.109) 

e^0 k ex e^0 k ex 

Then the trajectories satisfy the dynamical law 

in{t) = -2d i aV ± / 1 (a i (0) - 2d^ (Vh { ai {t)) - VVoM*))) • (4.110) 
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Proof. Fix a time t\ G (0, T*), an index i G {l,...,n}, and a constant < 5 < T* — tj. By 
the Holder continuity of tij and the energy bound ( I4.9p . we may find < T\ < r 2 < <7*/4 with 
r i = r i(£), r 2 = r 2(<5) so that cij(t) G S(aj(ti),ri) and 

/ ~9°e(ve,B £ )(t) < k ex C(5,n,Q) (4.111) 

for all t G + 5]. Now fix a unit vector e G 1R 2 and a function ^ G C%°(B(ai(ti), r 2 )) so that 

< i> < 1 and ^ = 1 on £(aj(ti), (ri+r 2 )/2). Define the function G C™(Q) by </>(a;) = (e-x)-0(x), 
and note that 0(aj(i)) = aj(t) ■ e and V0(aj(t)) = e for all t G + 5]- We may then apply 

Proposition 14. 141 to deduce that for any t 2 G \t\,t\ + S], 

vr(ai(t 2 ) - Oi(ti)) ■ e = -2^ fV Z ± {a i {t))dt + / / D 2 <j) : S(K), (4.112) 

where Z is defined in Proposition 14.141 It immediately follows that we can calculate the limit from 
the right: 

tt lim {ai(t) ~ a f l)} " 6 = -2nd ie ■ Z\ ai {h)) + [ D^: S{h',{ti)). (4.113) 

*->*! t ~~ *1 ■/B(oi(ti),ra)\B(o t (ti),n) 

A similar argument shows that the left limit exists and agrees with the right limit. 

In the annulus A := B(ai(ti), r 2 )\5(aj(ti), r\) the function h*(ti) is smooth and satisfies 
—Ah' m (ti) + K(ti) = 0. A direct calculation shows that for any smooth h, 

div 3(h) = (-Ah + h)Vh, (4.114) 

so divS'(^(ti)) = in A. Integrating by parts and using the structure of 0, this implies that 

D 2 <P : S(K(h)) = - [ e • S(K(h))u. (4.115) 

A J dB(pi(ti),ri) 

We may thus combine ( 14. 113ft (with t — >• t± replaced with the full limit t — > ti) and ( 14.1151) to 
deduce that a, is differentiable at t\ and that 

Tta^h) = -27rd i Z(o i (t 1 )) - / ^(^(ti))!/. (4.116) 

The second through fourth parameter regimes are easier to deal with since k ex 3> 1, which forces 
h'x = and S(h'^) = 0. To derive the dynamics, we must only determine the structure of Z 1 - in 
terms of the regime. It is easy to see that 

V" 1 /! in regime 2 

V/io-VVo in regime 3 (4.117) 

,aV ± /i + /3(V/i -V ± / ) in regime 4, 

from which (14.1101) follow immediately since t\ was arbitrary. For the regime j ex = h ex = 1 we note 
that by sending 5 — > we may also let r± — >■ 0. Applying Lemma [4.111 and taking the limit r± — > 
in ( 14. 1161) then yields 

= --V ai W d (a(ti)) - 2d i Z- L (a i (*i)), (4.118) 

7T 

with Z x =V i /i + V/i - V^/o, from which 04108]) follows. 

□ 
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A Well-posedness and regularity 

In this section we record some results on the well-posedness and a priori estimates for (11. 6p . 

Proposition A.l. The time-dependent Ginzburg-Landau equations (11 .6p are well-posed for all time, 
and the solutions are smooth. 

Proof. It is a simple matter to see that (u, A, $) solve (jl.6p if and only if (v, B, $) solve (I2.6p - fl2.8p . 
Using this reformulation of the problem, the problem is amenable to standard fixed-point techniques 
for solving semi-linear parabolic problems. A straightforward modification of the method employed 
in [HI] yields well-posedness. Smoothness follows from standard bootstrapping. □ 

We also record the following L°° bounds on u and V^w, which follow from a simple modification 
of Proposition 2.8 in [3T] . 

Lemma A. 2. Suppose that the initial data satisfy \\uq\\ Lx ^ < 1 and \\ ^ AqUoWl^^ < C/e. Then 

\\u(t)\\ Loc{n) <l for allt>0 (A.l) 

and 

l|VAu(f)|| L oo ( n) < C / e f° r al1 1 > 0- ( A - 2 ) 

B Static analysis of the Ginzburg-Landau energy 

In this appendix we will record some energy estimates for the static Ginzburg-Landau energy that 
are useful in the analysis of the dynamics. In particular, when we know the limit 

n 

\i e := /j,(u e ,A E ) = curl((m e , V A e u e ) + A E ) ->■ 27r^G^ ai , (B.l) 

i=i 

we will derive various estimates in terms of the vortex locations and degrees. Most of the results 
are variants of well-known ones, but cannot be found in the literature in the exact form we need. 
In addition to the free energy F £ we will also use the weighted free energy F£, which is given by 

F^u, *) = \J I Va«| 2 + ^(1 - M 2 ) 2 + r 2 |curl A\ 2 (B.2) 

for some r > 0. The simplified energy is given by E £ (u) = F e (u,0). We employ the notation 
F £ {u, A, S) for S C VL to mean the energy with the integral evaluated only over S. 

Our first result gives a lower bound for the free energy in terms of the limiting vorticity measure. 

Proposition B.l. Suppose (u E ,A e ) satisfy the bound F e {u £ ,A £ ) < C |loge| as well as ||w e || L «, < 1, 
IIVa^sII^oo < C/e, and Va s u £ ■ v — on dQ. Further suppose that 

n 

fi E ^27rJ2 d iK, (B.3) 

i=l 

where di = ±1, and that ||curl A £ || L2 < K £ for < K £ ^ |loge|. Then 

F £ {u £ , A £ ) > irn | log e | + ni + W d (a) + o(l)K £ . (B.4) 
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Proof. For convenience we will drop the subscript e in the proof, writing u, A, \i in place of u £ , A £ , \x e . 
Since we want a lower bound of F £ , we may first perform a minimization of F e (u,A) over all A, 
keeping u fixed. That is, we bound 

F e (u,A) > mm F £ (u, A) := F £ (u,B) (B.5) 

where B solves 

V" 1 curl 5 = — (iu, V bu) in £1 r> 

(B.o) 

curl B = on <9fi. 

We also fix the Coulomb gauge so that div B = in Q and B ■ v = on <9fi. By the Poincare lemma, 
this allows us to write B = V -1 ^, where £ = on 50. Since curll? = A£ and |Vcurl5| < |Vsf|, 
elliptic regularity gives that 

U\\ HH n) < \\cut\A\\ l2 <K £ , 

Mf H z(n)< c I \cm\A\ 2 + \V A v\ 2 <CF e {v,A) <C\\oge\. (B ' ?) 



We see from the first bound that ||£||c? > o (fi) — CK e f° r an y a ^ -0' anc ^ fr° m the second that 

^ < C | log e | . 

We expand F £ (u, B), employing the fact that B solves flB.6j) . to see that 



F £ (u, B) = E e (u) - - / \B\ 4 + \cuABy + o(l). (B.8) 



Taking the curl of ( 1B.6j) and adding curl B to both sides, we have that — A 2 £ + A£ = \i in Q, with 
£ = A£ = on <9fi. Then, integrating by parts, we have that 



\B\ 4 + Icurl^l" = - / |V£| 2 + |A£| 2 



^Jj(A^-AO (B.9) 



Define as the solution to 



4 / u>. 



-A 2 ^ + A^ = 2vrE^a l mfi 
£ = A& = on da 



Since £/K £ is bounded in C ' a (f2) for any a e (0, 1), we may assume that up to extraction 

jfe. ifif. = o(i) 

# £ 10 if l<iir e < |loge|. 



In either case 



„ n 

/ ^ = 2ttV^*K) + o(1)^- (B.12) 
Jn i=l 
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To deal with the E £ (u) term we have to show that the hypotheses carry over to u. First note 
that 

C C 

l|V«|| Loo < || V B u\\ Loo + \\u\\ Loa \\B\\ Loc <— + Cy/\loge\ < - (B.13) 

and that E e (v) < C |loge| . Since B ■ v = on the boundary, it holds that Vw • v = Vbv ■ v = on 
dQ. Finally, we note that since \x = curl((m, Vu) + (1 — \u\ 2 )B), the convergence of \x guarantees 
that curl(w, Vu) — > lix Y^h=i as wen - We ^ nus nave that all of the hypotheses of Proposition 
4.3 of [28] are satisfied; we find that 

E £ (u) > -rm [log e| + + — 7r djdj log |a» — %| + n'^2d i djRn(ai, dj) J + o(l), (B.14) 



where Rq is defined by Rn(x, y) = P(x, y) + log \x — y\ , with P is the solution to 

(B.15) 



- A x P(x, y) = 2Tc5 y in Q 



P(x, y) = for x € cftl 



We deduce that 
F £ (u, A) > im | log e | + wy 



-7r log |aj — a_j| + 7r ^iC*( a i) + /J didjRn(ai, a,j) J + o(l)lf e . (B.16) 

i^j i=l i,j / 

A simple calculation shows that 

n 

^ di^(ai) + ^2 didjR n (x, aA = ^ didjSa(ai, a 3 -). (B.17) 

i=l ij i,j 

Substituting this into (IB. 161) yields the desired lower bound. 

□ 

The next result gives a lower bound for the weighted free energy in balls near the vortex locations. 
Lemma B.2. Suppose that (u £ , A £ ) are such that 

Fe(u £ , A £ ) < im |loge| + o(|loge|) (B.18) 

and that 

n 

fi £ -> n = 27r^di5 ai , (B.19) 

i=l 

where di = ±1. Lei < a < | min{dist(aj, dQ)} U {|aj — Oj| | i ^ j}. Then there is a universal 
constant C > so that 

rm (log - - C) < F// 3 (u £ , A £ , U5(a i5 a)). (B.20) 
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Proof. We apply the ball construction (Theorem 4.1 of [27]) in each ball B^a^, a) to find a collection 
of balls Bi = {B it i, . . . , B itm -} with final radius r = |loge| -2 . For each i — 1, . . . , n and j = 1, . . . , m i; 
we let bij denote the center of the ball B^ and dB ij denote its degree. Then, according to the 



Theorem 6.1 of [27J, for each i = 1, . . . ,n it holds that 



E 

l<J<TOi 

BijCB(ai,a-e) 



dt 



< 



C 



\loge 



^F £ (u £ ,A £ 



oil). 



(B.21) 



(C°'\B( ai ,a))r 



From this and assumption flB.19j) . we deduce that for each i — 1, . . . , n there exists a ball A,j G Bi 
such that Bij C -B(aj, cr — e), gJb^- | > 1, and b^j — )■ a; as e — >• 0. By relabeling the collection Bi we 
may assume that Aj = A,i- 

The energy estimates of the ball construction imply that 



vrA lo. 



eD,, 



(B.22) 



where 



A 



E 



l<j<mi 
Bij<ZB(ai,cr-s) 



Moreover, Di < Cn for some universal constant C. Summing (IB.22[) over 



, n and plugging 



in the value of r and the upper bound flB.18j) . we find that 

ran. ra 

7T A | log e I - 2tt ^ Di log | log e | - tt A (log A + C) < nn | log £ | + o ( | log e | ) . ( B . 23 ) 



4=1 



Dividing by |loge|, we deduce that ^D^ <nfor e is sufficiently small. However, since for each 
i = l,...,n the ball A,i C B(ai,o~ — e) and |<i,B i:L | > 1, the reverse inequality n < ^2 Di must 
also hold. Hence E^ A = n. From this we deduce that A — 1 f° r each i and that . = for 
j = 2, . . . ,m r 

We now grow the balls in each collection Bi into a larger collection, B[, with total radius r = a/3. 
There must exist a ball B\ G B\ such that A,i C B[. By the above analysis on the degrees of the 
balls in Bi it holds that deg(A') = di and that the degrees of the other balls in B[ vanish. Moreover, 
since rad(A') < cr/3 and B\ 3 b^\ — > a i: it must be the case that B[ C B(ai, a — e). Plugging into 
the ball construction again then yields the bound 



7T 



(log |- - C) < F e CT / 3 (n e , A £ , B(a h a)) 



for each i = 1, . . . ,n, from which (1B.20|) immediately follows. 



(B.24) 



□ 



We now use the previous lemmas to prove that curl A £ is lower order than F £ and that the energy 
away from the vortex locations is bounded. 

Lemma B.3. Suppose (u £ ,A £ ) satisfy ji £ — > 2n Yji=i diO~ ai with di = ±1 as well as the bound 

F £ (u £ ,A £ )<irn\\oge\+K £ (B.25) 
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for some sequence < K £ <C | log e | . Let 



< a < - min{dist(a;, dQ)} U {|a 4 - aA \ i^ j} U {6/V2}. (B.26) 
Write Q a = Q\ U S(a i; a). Then 

If \V As u £ \ 2 + ^(l-\u £ \ 2 ) 2 + l [ \cuilA £ \ 2 <7mlog- + Cn + K £ , (B.27) 
•where Cq is a constant depending on Q and C is a universal constant. 

Proof. Write F £ (u £ ,A £ ) = F £ (u £ , A £ ,fl a ) + F £ (u £ , A £ ,UB(ai,a)). We further decompose the latter 
term to 

F £ {u £ ,A £ ,UB(a l ,a)) = F^ 3 (u £ ,A £ ,UB(a i ,a)) + ]-(l-^-) [ |curlA £ | 2 . (B.28) 

1 \ 9 / JuB(a t ,a) 

By construction 

\ (l - ^) [ |curM e | 2 >\l |curlA £ | 2 . (B.29) 

Lemma TB . 2 1 shows that F £ ^ 3 (u £ , A £ , UB(ai, a)) > nn (log - — C) for C a universal constant. Chain- 
ing together these lower bounds with the upper bound of the hypothesis then gives 

\f \V Ae u £ \ 2 + 7 ^(l-\u £ \ 2 ) 2 + ] [ \cm\A £ \ 2 <nnlog- + Cn + K E , (B.30) 
the desired inequality. □ 
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